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Abstract Involutive bases are a special form of non-reduced Grobner bases with 
additional combinatorial properties. Their origin lies in the Janet-Riquier theory of 
linear systems of partial differential equations. We study them for a rather general 
class of polynomial algebras including also non-commutative algebras like those 
generated by linear differential and difference operators or universal enveloping 
algebras of (finite-dimensional) Lie algebras. We review their basic properties us- 
ing the novel concept of a weak involutive basis and concrete algorithms for their 
construction. As new original results, we develop a theory for involutive bases 
with respect to semigroup orders (as they appear in local computations) and over 
coefficient rings, respectively. In both cases it turns out that generally only weak 
involutive bases exist. 



1 Introduction 

In the late 19th and early 20th century a number of French mathematicians de- 
veloped what is nowadays called the Janet-Riquier theory of differential equations 
||37ll38ll39ll48ll53]|57ll58l . It is a theory for general systems of differential equa- 
tions, i. e. also for under- and overdetermined systems, and provides in particular 
a concrete algorithm for the completion to a so-called passiv43 system. In recent 
times, interest in the theory has been rekindled mainly in the context of Lie symme- 
try analysis, so that a number of references to modern works and implementations 
are contained in the review ll36l . 

The defining property of passive systems is that they do not generate any non- 
trivial integrability conditions. As the precise definition of passivity requires the 
introduction of a ranking on the set of all derivatives and as every linear system of 

' Sometimes the equivalent term "involutive" is used which seems to go back to Lie. 
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partial differential equations with constant coefficients bijectively corresponds to a 
polynomial module, it appears natural to relate this theory to the algebraic theory 
of Grobner bases |[Tl|6l. 

Essentially, the Janet-Riquier theory in its original form lacks only the concept 
of reduction to a normal form; otherwise it contains all the ingredients of Grobner 
bases. Somewhat surprisingly, a rigorous links has been established only fairly 
recently first by Wu IMI and then by Gerdt and collaborators who introduced 
a special form of non-reduced Grobner bases for polynomial ideals f20 '21''62J, 
the involutive bases (Wu's "well-behaved bases" correspond to Thomas bases in 
the language of |20|). A slightly different approach to involutive bases has been 
proposed by Apel |4|; it will not be used here. 

The fundamental idea behind involutive bases (originating in the pioneering 
work of Janet I37II38I ) is to assign to each generator in a basis a subset of all vari- 
ables: its multiplicative variables. This assignment is called an involutive division, 
as it corresponds to a restriction of the usual divisibility relation of terms. We only 
permit to multiply each generator by polynomials in its multiplicative variables. 
As we will see later in Part II, for appropriately prepared bases, this approach 
automatically leads to combinatorial decompositions of polynomial modules. 

Like Grobner bases, involutive bases can be defined in many non-commutative 
algebras. We will work with a generalisation of the polynomial algebras of solvable 
type introduced by Kandry-Rodi and Weispfenning [41J. It is essentially equivalent 
to the generalisation discussed by Kredel 1421 or to the G-algebras considered by 
Apel [2 J and Levandovskyy ll44ll45ll . In contrast to some of these works, we explic- 
itly permit that the variables act on the coefficients, so that, say, linear differential 
operators with variable coefficients form a polynomial algebra of solvable type in 
our sense. Thus our framework automatically includes the work of Gerdt ifTSl on 
involutive bases for linear differential equations. 

This article is the first of two parts. It reviews the basic theory of involutive 
bases; this is immediately done in the framework of polynomial algebras of solv- 
able type, as it appears to be the most natural setting. Indeed, we would like to 
stress that in our opinion the core of the involutive bases theory is the monomial 
theory (in fact, we will formulate it in the language of multi indices or exponent 
vectors, i. e. in the Abelian monoid (Ng , +), in order to avoid problems with non- 
commuting variables) and the subsequent extension to polynomials requires only 
straightforward normal form considerations. 

While much of the presented material may already be found scattered in the 
literature (though not always in the generality presented here and sometimes with 
incorrect proofs), the article also contains some original material. Compared to 
Gerdt and Blinkov [20J , we give an alternative definition of involutive bases which 
naturally leads to the new notion of a weak involutive basis. While these weak 
bases are insufficient for the applications studied in Part II, they extend the appli- 
cability of the involutive completion algorithm to situations not covered before. 

The main emphasis in the hterature is on optimising the simple completion 
algorithm of Section|7]and on providing fast implementations; as the experiments 
reported in |23 1 demonstrate, the results have been striking. We will, however, ig- 
nore this rather technical topic and instead study in Part II a number of applications 
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of involutive bases (mainly for the special case of Pommaret bases) in the struc- 
ture analysis of polynomial modules. This will include in particular the relation 
between involutive bases and the above mentioned combinatorial decompositions. 
Note, however, that in these applications we will restrict to the ordinary commuta- 
tive polynomial ring. 

This first part is organised as follows. The next section defines involutive di- 
visions and bases within the Abelian monoid (Nq,+) of multi indices. It also 
introduces the two most important divisions named after Janet and Pommaret, re- 
spectively. Section [3] introduces the here used concept of polynomial algebras of 
solvable type. As the question whether Hilbert's Basis Theorem remains valid is 
non-trivial if the coefficients form only a ring and not a field, Section|4]is devoted 
to this problem. The following three sections define (weak) involutive bases and 
give concrete algorithms for their construction. 

The next four sections study some extensions of the basic theory. Section [8] 
analyses the relation between left and right ideals in polynomial algebras of solv- 
able type and the computation of bases for two-sided ideals; this requires only 
a straightforward adaption of classical Grobner basis theory. The following three 
sections contain original results. The first two ones generalise to semigroup orders 
and study the use of the Mora normal form. Finally, Section [TT| considers involu- 
tive bases over rings. It turns out that in these more general situations usually only 
weak bases exist. 

In a short appendix we fix our conventions for term orders which are inverse to 
the ones found in most textbooks on Grobner bases. We also prove an elementary 
property of the degree reverse lexicographic term order that makes it particularly 
natural for Pommaret bases. 



2 Involutive Divisions 

We study the Abelian monoid (Nq , +) with the addition defined componentwise 
and call its elements multi indices. They may be identified in a natural way with the 
vertices of an n-dimensional integer lattice, so that we can easily visualise subsets 
of Nq . For a multi index 1^ € Ng we introduce its cone C{i') = v + Ng , i. e. the 
set of all multi indices that can be reached from v by adding another multi index. 
We say that v divides /i, written v \ (x,\f ^ £ C{v). Given a finite subset N C Nq , 
we define its span as the monoid ideal generated by M: 



The basic idea of an involutive division is to introduce a restriction of the cone 
of a multi index, the involutive cone: it is only allowed to add multi indices certain 
entries of which vanish. This is equivalent to a restriction of the above defined 
divisibility relation. The final goal will be having a disjoint union in ^ by using 
only these involutive cones on the right hand side. This will naturally lead to the 
combinatorial decompositions discussed in Part II. 




(1) 
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In order to finally give the definition of an involutive division, we need one 
more notation: let iV C {1, . . . , n} be an arbitrary subset of the set of the first n 
integers; then we write = {i^ e Nq \ yj ^ N : i/j = O] for the set of all 
multi indices where the only entries who may be non-zero are those whose indices 
are contained in N. 

Definition 2.1 An involutive division L is defined on the Abelian monoid (Nq , +), 
if for any finite set M C Nq a subset Nl,j^{v) C {1, . . . , n} of multiplicative 
indices is associated to every multi index v ^ M such that the following two 
conditions on the involutive cones ClmW) + N^^ ^j.^^ are satisfied. 

1. If there exist two elements ii,!/ E M with Ci^7v(/i) H Ci jv(z^) ^ 0, either 
CL,j\r{tJ') C CL,Af{'^) orCL,Af{i^) ^ Clx{p) holds. 

2. IfM' C TV, then Nl,x{v) C NL^j^,{v)for all v £ TV'. 

An arbitrary multi index /i G Nq is involutively divisible by v £ M, written 

Before we discuss the precise meaning of this definition and in particular of 
the two conditions contained in it, we should stress the following important point: 
as indicated by the notation, involutive divisibility is always defined with respect 
to both an involutive division L and a fixed finite set TV C Nq : only an element of 
TV can be an involutive divisor. Obviously, involutive divisibility |l,aA M implies 
ordinary divisibility v \ [i. 

The involutive cone CL,is!iy) of any multi index e TV is a subset of the 
full cone C(y). We are not allowed to add arbitrary multi indices to v but may 
increase only certain entries of v determined by the multiplicative indices. The 
first condition in the above definition says that involutive cones can intersect only 
trivially: if two intersect, one must be a subset of the other. 

The non-multiplicative indices form the complement of Nl,^{v) in {1, . . . , n} 
and are denoted by Nl_j\/{v). If we remove some elements from the set TV and 
determine the multiplicative indices of the remaining elements with respect to the 
subset TV', we obtain in general a different result than before. The second condition 
for an involutive division says that while it may happen that a non-multiplicative 
index becomes multiplicative for some v G TV', the converse cannot happen. 

Example 2.2 A classical involutive division is the Janet division J. In order to 
define it, we must introduce certain subsets of the given set TV C Nq : 

(dfe, . . . ,dn) = [v £ Af \ Vi = di , k < i < n] . (2) 

The index n is multiplicative for v G TV, if i/„ = maxy^gyv" {Mk}, and fc < ti is 
multiplicative for v G (dfe+i, . . . , d„), if Vk = va&yi^^(^dk+i,...,d„) {A^fc}. 

Obviously, this definition depends on the ordering of the variables xi, . . . , a;„ 
and we obtain variants by applying an arbitrary but fixed permutation tt G S'„ to 
the variables. In fact, Gerdt and Blinkov [20] use an "inverse" definition, i. e. they 

first apply the permutation ^ '^^ ■ convention is the original one of 

Janet [39. pp. 16-17]. 
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Gerdt et al. Il22l designed a special data structure, the Janet tree, for the fast 
determination of Janet multiplicative indices and for a number of other operations 
useful in the construction of Janet bases (Blinkov [lOj discusses similar tree struc- 
tures also for other divisions). As shown in ll32ll . this data structure is based on a 
special relation between the Janet division and the lexicographic term order (see 
the appendix for our non-standard conventions). This relation allows us to com- 
pute very quickly the multiplicative variables of any set M with Algorithm[T] The 
algorithm simply runs two pointers over the lexicographically ordered set M and 
changes accordingly the set M of potential multiplicative indices. <\ 



Algorithm 1 Multiplicative variables for the Janet division 

Require: finite list A/" = {z^'^^ . . . , i^''^' } of pairwise different multi indices from 
Ensure: list N = {Nj^u{v^^^), . . . , Nj^j^f{u''''^)} of fists with multiplicative variables 

V. Af ^ sort(A/', ^lex); I' ^ ^f['^] 

2: pi^n; M <- {1, . . . , n}; N[l] «- M 

3: for j from 2 to \Af\ do 

4: p2^ma.x{i\{u-J\f[j])^^0}- M^M\{p2} 

5: if pi < p2 tiien 

6: M ^ MLl{pi,...,p2-l} 

7: end if 

8: N[j]^M; u ^ f^[j]; pi ^ P2 
9: end for 
10: return iV 



Definition 2.3 The division L is globally defined, ;/ the assignment of the multi- 
plicative indices is independent of the set J\f; in this case we write simply Nl{v). 

Example 2.4 Another very important division is the Pommare^ division P. It as- 
signs the multiplicative indices according to a simple rule: if 1 < fc < n is the 
smallest index such that Vk > for some multi index i/ G Nq \ {[0, . . . , 0]}, then 
we call k the class of written els z/, and set Np{y) = {1, . . . , fc}. Finally, we 
define iVp([0, . . . , 0]) = {1, . . . , n}. Hence P is globally defined. Like the Janet 
division it depends on the ordering of the variables xi, . . . , a;„ and thus one may 
again introduce simple variants by applying a permutation. 

Above we have seen that the Janet division is in a certain sense related to the 
inverse lexicographic order. The Pommaret division has a special relation to class 
respecting orders (recall that according to Lemma PV. II anv class respecting term 
order coincides on terms of the same degree with the reverse lexicographic order). 
Obviously, for homogeneous polynomials such orders always lead to maximal sets 
of multiplicative indices and thus to smaller bases. But we will also see in Part II 
that from a theoretical point of view Pommaret bases with respect to such an order 
are particularly useful. < 

^ Historically seen, the terminology "Pommaret division" is a misnomer, as this division 
was already introduced by Janet |37, p. 30], too. However, the name has been generally 
accepted by now, so we stick to it. 
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Above we introduced the span of a set A/^ C Ng as the union of the cones of its 
elements. Given an involutive division it appears natural to consider also the union 
of the involutive cones. Obviously, this yields in general only a subset (without 
any algebraic structure) of the monoid ideal (Af). 



Definition 2.5 The involutive span of a finite set Af C Ng is 

{M)l = U CmA^) ■ 



(3) 



The set N is weakly involutive /or the division L or a weak involutive basis of the 
monoid ideal {M), if {M) l = (-A/"). A weak involutive basis is a strong involutive 
basis or for short an involutive basis, if the union on the right hand side of ^ is 
disjoint, i. e. the intersections of the involutive cones are empty. We call any finite 
set Af N C Ng such that (Af) l ~ {Af) a (weak) involutive completion of Af. 
An obstruction to involution /or the set Af is a multi index v G (Af) \ {Af)h- 

Remark 2.6 An obvious necessary condition for a strong involutive basis is that no 
distinct multi indices fi^v E N exist such that /i \ v. Sets with this property are 
called involutively autoreduced. One easily checks that the definition of the Janet 
division implies that C^fj{^) D Cjsf^L{v) = whenever [i ^ v. Hence for this 
particular division any set is involutively autoreduced. <l 



• 
• 


• • • 


<l 


• • • 



[2,0] 




Fig. 1 Left: intersecting cones. Right: involutive cones. 



Example 2. 7 Figure [T] demonstrates the geometric interpretation of involutive di- 
visions for n = 2. In both diagrams one can see the monoid ideal generated by the 
set Af = { [0, 2] , [2, 0] } ; the vertices belonging to it are marked by dark points. The 
arrows represent the multiplicative indices, i. e. the "allowed directions", for both 
the Janet and the Pommaret division, as they coincide for this example. The left 
diagram shows that the full cones of the two elements of Af intersect in the darkly 
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shaded area and that M is not (weakly) involutive, as the muhi indices [fc, 1] with 
A: > 2 are obstructions to involution. The right diagram shows a strong involutive 
basis of {M) for both the Janet and the Pommaret division. We must add to M the 
multi index [2, 1] and both for it and for [2, 0] only the index 1 is multiplicative. 
One clearly sees how the span {M) is decomposed into three disjoint involutive 
cones: one of dimension 2, two of dimension 1. <\ 

We are particularly interested in strong involutive bases. The following result 
shows that in the "monomial" case any weak involutive basis can be reduced to a 
strong one by simply eliminating some elements. 

Proposition 2.8 If M is a weakly involutive set, then a subset J\f' C AA exists such 
that M' is a strong involutive basis of (A/"). 

Proof This proposition represents a nice motivation for the two conditions in Def- 
inition |2T| of an involutive division. If M is not yet a strong involutive basis, the 
union in Q is not disjoint and intersecting involutive cones exist. By the first con- 
dition, this implies that some cones are contained in other ones; no other form of 
intersection is possible. If we eliminate the tips of these cones from JV, we get a 
subset AA' C M which, by the second condition, has the same involutive span, as 
the remaining elements may only gain additional multiplicative indices. Thus after 
a finite number of such eliminations we arrive at a strong involutive basis. □ 

Remark 2.9 Let Zi, X2 be two monoid ideals in Nq and TVi, A/2 (weak) involutive 
bases of them for some division L. In general, we cannot expect that A/i U A/2 is 
again a weak involutive basis of the ideal Zi + Z2, as the involutive cones of the 
generators may shrink when taken with respect to the larger set Afi U A/2. Only 
for a global division we always obtain at least a weak involutive basis (which may 
then be reduced to a strong basis according to Proposition l2.8b . < 

Recall that for arbitrary monoid ideals a basis JV is called minimal, if it is not 
possible to remove an element of Af without losing the property that we have a 
basis. A similar notion can be naturally introduced for involutive bases. 

Definition 2.10 Let X C Nq be a monoid ideal and L an involutive division. An 
involutive basis Af of I with respect to L is called minimal, if any other involutive 
basis Af' of X with respect to L satisfies Af C TV'. 

Obviously, the minimal involutive basis of a monoid ideal is unique, if it exists. 
For globally defined divisions, any involutive basis is unique. 

Proposition 2.11 Let L be a globally defined division and X C Nq a monoid 
ideal. If X has an involutive basis for L, then it is unique and thus minimal. 

Proof Let Af be the minimal basis of X and Afi, N2 two distinct involutive bases 
of X. Both Afi \ N2 and Af2 \ Afi must be non-empty, as otherwise one basis 
was contained in the other one and thus the larger basis could not be involutively 
autoreduced with respect to the global division L. Take an arbitrary multi index 
v G A/i \ A/2. The basis Af2 contains a unique multi index /x such that /i |i, i^. 
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It cannot be an element of A^i, as A/i is involutively autoreduced. Thus A/i must 
contain a unique multi index A 7^ /i such that X\l fi. As L is globally defined, this 
implies that A | ^ and A ^ i', a contradiction. □ 

The algorithmic construction of (weak) involutive completions for a given set 
M C INg will be discussed in detail in Section|6l For the moment we only note that 
we cannot expect that for an arbitrary set J\f and an arbitrary involutive division L 
an involutive basis AA' of (JV) exists. 

Example 2.12 We consider the set A/" = {[1,1]} for the Pommaret division. As 
els [1, 1] = 1, we get Np[[l, 1]) = {1}. So Cp([l, 1]) C C([l, 1]). But any multi 
index contained in {M) also has class 1. Hence no finite involutive basis of {M) 
exists for the Pommaret division. We can generate it involutively only with the 
infinite set {[1, fc] | fc e N}. 

More generally, we may consider an irreducible monoid ideal I in Ng . It is 
well-known that any such I has a minimal basis of the form {(^i)ij^ , . . . , {ik)ik } 
with \ < k < n, £j > Q and 1 < ii < ■ ■ ■ < ik < n. Here is the multi index 
where all entries are zero except of the ijth one which has the value Ij. Such an 
ideal possesses a Pommaret basis, if and only if there are no "gaps" in the sequence 
ii < • ■ • < ifc < n, i. e. ifc = n and ii = n — fc + 1. Indeed, if a gap exists, say 
between ij and ij+i, then any Pommaret basis must contain the infinitely many 
multi indices of the form + (^)i +1 with ^ > and thus cannot be finite 

(obviously, in this case a simple renumbering of the variables suffices to remedy 
this problem). Conversely, if no gaps appear, then it is easy to see that the set of all 
multi indices [0, . . . , 0, li. , Hi.+i, . . . , /i„] with 1 < j < fc and < /i^ < in-k+i 
is a strong Pommaret basis of X. < 

Definition 2.13 An involutive division L is called Noetherian, if any finite subset 
J\f C Ng possesses a finite involutive completion with respect to L. 

Lemma 2.14 The Janet division is Noetherian. 

Proof Let J\f C Ng be an arbitrary finite set. We explicitly construct a Janet basis 
for (AA). Define the multi index /i = 1cm AA by /i; = max^g^ Vi. Then we claim 
that the (obviously finite) set 

AT- {PG (AA) l^eC(P)} (4) 

is an involutive completion of AA with respect to the Janet division. Indeed, AA C AA 
and AA C (TV) . Let p G (AA) be an arbitrary element. If p G AA, then trivially 
p G (-AA) J. Otherwise set I — {i \ Pi > fJ,i} and define the multi index p by 
Pi — Pi for i ^ I and pi = pi for i G I, i. e. pi = min {pi, /ij}. By construction 
of the set AA and the definition of p, we have that p G M and / C Nj fj{p). But 
this implies that p G C j j^{p) and thus AA is a finite Janet basis for (AA). □ 
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3 Polynomial Algebras of Solvable Type 

We could identify multi indices and monomials and proceed to define involutive 
bases for polynomial ideals. But as the basic ideas remain unchanged in many dif- 
ferent situations, e. g. rings of linear differential or difference operators, we gen- 
eralise a concept originally introduced by Kandry-Rody and Weispfenning BTl 
and use polynomial algebras of solvable type (the same kind of generalisation has 
already been intensively studied by Kredel (42]). 

Let V = TZ[xi, . . . , Xn] be a polynomial ring over a unitary rin^ TZ. If TZ 
is commutative, then T' is a unitary commutative ring with respect to the usual 
multiplication. We equip the 7?.-module V with alternative multiplications, in par- 
ticular with non-commutative ones. We allow that both the variables Xi do not 
commute any more and that they operate on the coefficients. The usual multiplica- 
tion is denoted either by a dot • or by no symbol at all. Alternative multiplications 
V X V ^ V are always written as f * g. 

Like Grobner bases, involutive bases are always defined with respect to a term 
order -<. It selects in each polynomial / G "P a leading term f — x^ with 
leading exponent le^ f = ji. The coefficient r G 7^ of a;^ in / is the leading 
coefficient \c-^ f and the product rx^ is the leading monomial Im^ /. Based on 
the leading exponents we associate to each finite set c "P a set le^ T C Nq to 
which we may apply the theory developed in the previous section. But this requires 
a kind of compatibility between the multiplication -k and the chosen term order. 

Definition 3.1 {V, *, ■<) is a polynomial algebra of solvable type over the unitary 
coefficient ring TZ for the term order -<, if the multiplication -k : V x V ^ V 
satisfies three axioms. 

(i) {P, k) is a ring with unit 1. 

(ii) yr€n,f£P: r*f = rf. 

(iii) V/x,:/ G Nff, r G 7^\ {0} : le^ {xt" * x") = fi + ly A le^ (a;^*r) = ^. 

Condition (i) ensures that arithmetics in (P,*, -<) obeys the usual associative 
and distributive laws. Because of Condition (ii), {V, *, -<) is a left 7?.-module. We 
do not require that it is a right 7?,-module, as this would exclude the possibility 
that the variables Xi operate non-linearly on TZ. Condition (iii) ensures the com- 
patibility of the new multiplication ★ and the term order we say that the mul- 
tiplication ★ respects the term order ^. It implies the existence of injective maps 
Pp, : TZ ^ TZ, maps -.TZ ^ T' with le^ (ft,^(r)) ^ p. for all r £ TZ, coefficients 
ffiiy G 7?. \ {0} and polynomials /i^i, G P with le^ /i^^ ^ fj, + u such that 

xf'-kr:^ p^{r)x'' + h^{r) , (5a) 
x^ * x" ^ r^,x''+'' + h^, . (5b) 

Lemma 3.2 The maps and the coefficients r^^ satisfy for arbitrary multi in- 
dices Hjiy, X € Ng and for arbitrary ring elements r TZ 

Pt^{Piy{r))rf„y = rf,t,Pf,+t,{r) , (6a) 
Pti{r^x)rf_i,u+x = rf_,urfj,+u,x . (6b) 



^ For us a ring is always associative. 
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Furthermore, all maps are ring endomorphisms. 

Proof The first assertion is a trivial consequence of the associativity of the muhi- 
plication -k. The equations correspond to the leading coefficients of the equalities 
a;^ * {x'^ * r) = (x^ ★ x") * r and * {x'^ ★ x^) — (x^ -kx") -k x^, respectively. 
The second assertion follows mainly from Condition (i). □ 

If 7?. is a (skew) field, then for arbitrary polynomials f,g G 7-* an element 
r e 7?. \ {0} and a polynomial h ^ V satisfying le^ h -< le^ (/ • g) exist such that 

f^g^rif ■g) + h. (7) 

Indeed, if Im^ / — ax^ and Im^ g = hx'^ , then a simple computation yields that 
r is the (unique) solution of the equation ap^{b)r^i, ~ rah and h is the difference 
f -k g — r{f ■ g). Under this assumption we may reformulate (iii) as 

(iii)' V/, g^V: le^ (/ * g) = le^ / + le^ g. 

For this special case the same class of non-commutative algebras was introduced in 
012II13I under the name PBW algebras (see Example l3 .5 I below for an explanation 
of this name). 

Proposition 3.3 The product -k is fixed, as soon as the following data are given: 
constants Vij G TZ\{0}, polynomials hij G V and maps pi : TZ ^ TZ, hi : TZ ^ V 
such that for 1 < i < n 

Xi*r ^ Pi{r)xi + hi{r) , Vr G 7^ , (8a) 

Xi -k Xj — TijXj -kXi + hij , VI < j < z . (8b) 

Proof The set of all "monomials" x^^ ★ x^^ ★ • • • ★ Xi with ii < 12 < ■ • • < 
forms an 7?.-linear basis of V, as because of (iii) the 7?.-linear map defined by 
Xi^ -k Xi^ -k ■ ■ ■ -k Xi ^ Xij ■ Xi^ - ■ ■ Xi is an 72.-module automorphism mapping 
the new basis into the standard basis. Obviously, it is possible to evaluate any 
product / *5 by repeated applications of the rewrite rules (O provided / and g are 
expressed in the new basis. □ 

Note that this proof is non-constructive in the sense that we are not able to 
determine the multiplication in terms of the standard basis, as we do not know ex- 
plicitly the transformation between the new and the standard basis. The advantage 
of this proof is that it is valid for arbitrary coefficient rings TZ. Making some as- 
sumptions on TZ (the simplest possibility is to require that it is a field), one could 
use Lemma[32]to express the coefficients r^^, and in (|5]l by the data in (O. This 
would yield a constructive proof. 

These considerations also show the main difference between our definition 
of solvable algebras and related definitions that one can find at various places in 
the literature. Our Condition (iii) represents in most approaches a lemma. Instead 
one usually imposes conditions like for i < j. Then the rewrite 

rule (l8b] l suffices to obtain an explicit transformation between the two bases used 
in the proof above. This is for instance the approach taken by Kandry-Rody and 
Weispfenning 11411 and subsequently by Kredel ll42l . 
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Of course, the data in Proposition 13.31 cannot be chosen arbitrarily. Besides 
the obvious conditions on the leading exponents of the polynomials hij and hi (r) 
imposed by Condition (iii), each map pi must be an injective 7^-endomorphism 
and each map hi must satisfy hi{r + s) = hi{r) + hi{s) and a kind of pseudo- 
Leibniz rule hiirs) ~ Pi{r)hi{s) + hi{r) ★ s. The associativity of ★ imposes 
further rather complicated conditions on the data. For the case of a G-algebra with 
the multiplication defined by rewrite rules they have been explicitly determined by 
Levandovskyy r44"45l who called them non-degeneracy conditions (see also the 
extensive discussion in [42, Sect. 3.3]). 

Example 3.4 An important class of non-commutative polynomials was originally 
introduced by Noether and Schmeidler |50| and later systematically studied by 
Ore ll52l : our exposition follows ifTTI . It includes in particular linear differential 
and difference operators (with variable coefficients). 

Let F be an arbitrary commutative ring and cr : F ^ F an injective endo- 
morphism. A pseudo-derivation with respect to a is a map (5 : F — > F such that 
(i) S{f + g) = Sif) + 6ig) and (ii) 6if ■g)=a{f)- S{g) + S{f)-g for all /, 5 £ F. 
If (T = idp, the identity map, (ii) is the standard Leibniz rule for derivations. If 
a ^ id-F, one can show that there exists an /i e F such that S — h{(j — idp). And 
conversely, if 5 ^ 0, there exists an /i e F such that a = hS + id]p. Ore called 
(t(/) the conjugate and 6{f) the derivative of /. 

Given a and S, the ring F[9;(t, (5] of univariate Ore polynomials consists of 
all formal polynomials in d with coefficients in F, i. e. of expressions of the form 
9 ~ ^21=0 -fi^^ wilh fi E F and q G INq. The addition is defined as usual. The 
variable d operates on an element / G F according to the rule 

d*f^<j{f)d + S{f) (9) 

which is extended associatively and distributively to define the multiplication in 
F[(?; cr, 6]: given two elements 61,62 G F[(9; a, 6], we can transform the product 
61 -k 62 to the above normal form by repeatedly applying (|9]l. The injectivity of the 
endomorphism cr ensures that deg {61 * ^2) = deg 61 + deg 62- We call F[9; a, 6] 
the Ore extension of F generated by a and S. 

A simple concrete example is given by choosing for F some ring of differ- 
entiable functions in the real variable x, say 

yielding linear ordinary differential operators with polynomial functions as co- 
efficients (i. e. the Weyl algebra over Q). Similarly, we obtain linear recurrence 
and difference operators. We set F = <C{n), the space of sequences with com- 
plex elements, and take for a the shift operator, i. e. the automorphism mapping s„ 
to Sn+i. Then A = a — idp is a pseudo-derivation. ]F[E;a,0] consists of linear 
ordinary recurrence operators, W[E; a, A] of linear ordinary difference operators. 

For multivariate Ore polynomials we take a set = {di, cr„} of F- 
endomorphisms and a set Z\ = {61, . . . ,Sn} where each Si is a pseudo-derivation 
with respect to (7i. For each pair {at, Si) we introduce a variable di satisfying a 
commutation rule (|9]l. If we require that all the maps ctj, Sj commute with each 
other, i. e. ai o Uj = aj o o-j. Si o Sj — Sj o Si and ai o Sj — Sj o ai for all i ^ j, 
one easily checks that di ★ dj = dj * di, i.e. the variables di commute. Setting 
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V — {9i , . . . , we denote hyW[D]S,A] the ring of multivariate Ore polyno- 
mials. Because of the commutativity of the di we may write the terms as 9^ with 
multi indices /x e Nq , so that it indeed makes sense to speak of a polynomial ring. 
Comparing with Proposition 13. 31 we see that we are in the special case where the 
maps hi always yield constant polynomials and the variables Xi commute. 

Finally, we show that (F[P; S, A] , is an algebra of solvable type for any 
term order -<. The product of two monomial operators and hd'^ is given by 

ad''irbd''= (^\acr^{S''{b))d^+'' (10) 

where (^) is a shorthand for Y[^=i (aO' '^^ ^ ''^i ^ o • • • o cr^" and similarly for 6. 
By the properties of a term order this implies 

le^ {ad" ★ bd") = /i + = Ic^ (ad") + le^ (bd") , (11) 

as any term d^^"^ appearing on the right hand side of ( fTOb divides 9''+'^ and thus 
^A+i/ ^ gp+i/ fQj. ^jjy jgj-jj^ order ^. <1 

Example 3.5 Bell and Goodearl Q introduced the Poincare-Birkhoff-Witt exten- 
sion (for short PBW extension) of a ring 7?. as a ring V ^ TZ containing a finite 
number of elements xi, . . . ,Xn G V such that (i) "P is freely generated as a left 
7?.-module by the monomials x" with /i G Nq, (ii) xi -k r — r Xi e 7^ for all 
r £ TZ and (iii) £71 + TZxi + ■ ■ ■ TZxn- Obviously, any such 

extension is a polynomial algebra of solvable type in the sense of Definition 13.11 
for any degree compatible term order. Other term orders generally do not respect 
the multiplication in V. 

The classical example of such a PBW extension is the universal enveloping 
algebra of a finite-dimensional Lie algebra g which also explains the name: 
the Poincare-Birkhoff-Witt theorem asserts that the monomials form a basis of 
these algebras i60l . They still fit into the framework developed by Kandry-Rody 
and Weispfenning PTI . as the Xi do not act on the coefficients. This is no longer 
the case for the more general skew enveloping algebras 9^#ii(g) where is a 
k-algebra on which the elements of g act as derivations ll47l Sect. 1.7.10]. <1 

Example 3.6 In all these examples, the coefficients j-^^ appearing in (|5]l are one; 
thus dSbl i are classical commutation relations. This is no longer true in the quan- 
tised enveloping algebras il/i(g) introduced by Drinfeld 1 17] and Jimbo pOl- For 
these algebra it is non-trivial that a Poincare-Birkhoff-Witt theorem holds; it was 
shown for general Lie algebras g by Lusztig [46]. Berger [8] generalised this re- 
sult later to a larger class of associative algebras, the so-called q-algebras. They are 
characterised by the fact that the polynomials hij in (l8b] i are at most quadratic with 
the additional restriction that hij may contain only those quadratic terms XkXi that 
satisfy i < k < £ < j and k — i = j — £. Thus any such algebra is a polynomial 
algebra of solvable type for any degree compatible term order. 

A simple concrete example is the q-Heisenberg algebra for a real q > (and 
q ^ 1). Let / be a function of a real variable x lying in some appropriate function 
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space. Then we introduce the operators 




Tqfix) = f{qx) , Xf{x) 



xf{x) . 



(12) 



It is straightforward to verify that these three operators satisfy the following q- 
deformed form of the classical Heisenberg commutation rules 

Sq-k X — X -k 5q + Tq , Sq -k Tq — qTq * (5g , Tq -k X = qX -k Tq . (13) 

Hence the algebra 'k[5q^Tq,x] is a polynomial algebra of solvable type for any 
degree compatible term order (but also for any lexicographic order with Tq -< 5q 
and Tq -< x). <\ 

Example 3.7 Let (P, be a polynomial algebra of solvable type with a degree 
compatible term order Then P is a filtered ring with respect to the standard 
filtration Sq = ® Vi and we may introduce the associated graded algebra by 
setting {gr^V)q = Uq/Uq-i. It is easy to see that gr^V is again a polynomial 
algebra of solvable type for ^. If in (O deg /i,;(r) = 0, dcg hij < I, pi ~ id-jz and 
rij = 1 (which is for example the case for all Poincare-Birkhoff-Witt extensions), 
then in fact gr-^V — (V,-), the commutative polynomial ring. Such algebras are 
sometimes called almost commutative ||47l Sect. 8.4.2]. < 

Proposition 3.8 // the ring TZ is an integral domain, then any polynomial algebra 
(■p, of solvable type over it is an integral domain, too, and a left Ore domain. 

Proof The first assertion is a trivial consequence of (Q: if Tl has no zero divisors, 
then f ■ g implies / ★ g ^ 0. Hence V does not contain any zero divisors. 

For the second one we must verify the left Ore conditions II15II51I : we must 
show that one can find for any two polynomials /, g g V with f * g ^ two 
further polynomials G V \ {0} such that (p k f ~ ip k: g. We describe now a 
concrete algorithm for this task. 

We set = {f,g} and choose coefficients j-q, sq G TZ such that in the dif- 
ference rog k f — sof -k g = hi the leading terms cancel. Then we perform a 
pseudoreduction of hi with respect to JTq. It leads with an appropriately chosen 
coefficient to G TZto an equation of the form 



where the remainder hi satisfies le^ hi ^ (le^ J-q). If hi — 0, we are done and the 
polynomials cf) = to''o5 ~ 4'o and ip — tf)SQf + ipQ form a solution of our problem. 
By Part (iii) of Definition 13. II we have le^ hi ~< le^ / + le^ g. This implies by 
the monotonicity of term orders that le^ 0o le^ g and le^ ipo ~< le^ /. Thus we 
have found a non-trivial solution. 

Otherwise we set J-i = L) {hi} and choose coefficients ri,si G TZ such 
that in the difference rif -k hi — sihi * / = /12 the leading terms cancel. Now 
we perform a pseudoreduction of /12 with respect to Ti. This computation yields 
a coefficient ti G TZ and polynomials , i/ji , pi G V such that 



tohi = 00 * / + V'o * 3 + ^1 



(14) 



tih2 ^ (pi k f + ipi -k g + pi -k hi + h2 



(15) 
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where the remainder /12 satisfies le^ ^ (le^ Ti).lf h2 = 0, then we are done, 
as we can substitute hi from (fT4l i and obtain thus for our problem the solution cf) — 
{hrif - pi)-^{hrQg~(f)Q)-tiSihi + (t)i and?/; = (iiri/-pi)*(toSo/-^o)+^i- 
By the same reasoning on the leading exponents as above, it is a non-trivial one. 

Otherwise we iterate: we set IF2 ^ {^2], choose coefficients r2, S2 
such that in the difference 7-2/ * /i2 — S2/12 * / = ^3 the leading terms cancel, 
compute the remainder /13 of a pseudo reduction of with respect to J^2 and 
so on. If the iteration stops, i.e. if the remainder vanishes for some value 
iV £ N, then we can construct non-zero polynomials (f), -0 with / = ^-kg hy 
substituting all remainders hi by their defining equations. The iteration terminates 
by a simple Noetherian argument: (le^ Tq) C (le^ Ti) C (le^ T2) C • • • is a 
strictly ascending chain of monoid ideals in Ng and thus cannot be infinite. □ 

Obviously, we can show by the same argument that 'P is a right Ore domain. 
The Ore multipliers 0, t/j constructed in the proof above are not unique. Instead of 
always analysing differences of the form rif -khi — Sihi * /, we could have used 
differences of the form rig -k hi — Sihi * g 01 we could have alternated between 
using / and g and so on. In general, each ansatz will lead to different multipliers. 

We have given here a direct and in particular constructive proof that V satisfies 
the left and right Ore conditions. Instead we could have tried to invoke Theo- 
rem 2.1.15 of BtI stating that any right Noetherian integral domain is also a right 
Ore domain. However, as we will see in the next section, if the coefficient ring TZ 
of V is not a field, then the question whether or not V is (left or right) Noetherian 
becomes nontrivial in general. 

Example 3.9 In the commutative polynomial ring one has always the trivial so- 
lution ip — g and ip = f. One might expect that in the non-commutative case 
one only has to add some lower terms to it. However, this is not the case. Con- 
sider the universal enveloping algebra of the Lie algebra so(3). We may write it as 
il(so(3)) — h[xi, X2, X3] with the multiplication* defined by the relations: 

Xl-kX2 = X1X2 , X2-*cXi — X1X2 — X3 , 

xi ★ = xixa , xs-k xi = xix^ + X2 , (16) 

X2-kX3 = X2X3 , X3*X2 = X2X3 - Xi . 

This multiplication obviously respects any degree compatible term order but not 
the lexicographic order. Choosing f = xi and g = X2, possible solutions for 
(p -k f = Ip -k g are ip — X2 — 1 and ip = X1X2 — 2x3 01 (p — X1X2 + X3 and 
ip = Xi — 1. They are easily constructed using the algorithm of the proof of 
Proposition 13.81 once with / and once with g. Here we must use polynomials of 
degree 2; it is not possible to find a solution of degree 1. <] 

4 Hilbert's Basis Theorem for Solvable Algebras 

A classical property of the ordinary polynomial ring V ~ TZixi, . . . , Xn], which is 
crucial in the theory of Grobner bases, is Hilbert's Basis Theorem. For our more 
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general class of polynomial algebras, it remains true only under additional assump- 
tions. As V is generally non-commutative, we must distinguish left, right and two- 
sided ideals and thus also study separately whether V is left or right Noetherian. 

With the exception of Section [8] we will exclusively work with left ideals and 
thus do not introduce special notations. This restriction to left ideals is not only for 
convenience but stems from the fundamental left-right asymmetry of Definition l3.1l 
of a polynomial algebra of solvable type where products r*a;'' and x''*r are treated 
completely different. For this reason we discuss only the question when V is left 
Noetherian (see also Remark l4~8] below). 

Most classical proofs of Hilbert's Basis Theorem consider only the univariate 
case and then extend inductively to an arbitrary (but finite) number of variables. 
However, this inductive approach is not possible in arbitrary polynomial algebras 
of solvable type, as the multiplication * does not necessarily restrict to a subal- 
gebra with fewer variables. A simple counterexample is provided by the universal 
enveloping algebra il(5o(3)) introduced in Example 13.91 where * cannot be re- 
stricted to the subspace ]k;[xi,a;2] since X2 * xi — xiX2 — X3. This observation 
motivates the following definition. 

Definition 4.1 (P,*, -<) is an iterated polynomial algebra of solvable type, if (i) 
V — Tl.[xi\[x2] ■ ■ ■ [xn] where each intermediate ring Vk = 'R-[xi\[x2] ■ ■ ■ [xk\ is 
again solvable for the corresponding restrictions of the multiplication ★ and the 
term order -< and(ii) we have the equality Xk^Vk-i+Vk-i = Vk-i^Xk + Vk-i- 

For iterated polynomial algebras of solvable type we may apply the usual in- 
ductive technique for proving a basis theorem, but we must still impose further 
conditions on the multiplication *. The following result is proven in |47 The- 
orem 1.2.9] for Ore algebras, but it is a trivial exercise to verify that the proof 
remains valid for more general algebras. 

Theorem 4.2 If (V, ~<) is an iterated polynomial algebra of solvable type over 
a left Noetherian ring TZ, then V is a left Noetherian ring, too. 

The second condition in Definition 14. 1 1 cannot be omitted, if a basis theorem 
is to hold. McConnell and Robson |47 Example 1.2. 11] provide a concrete coun- 
terexample of a univariate polynomial ring of solvable type which violates this 
condition and which is neither left nor right Noetherian. 

With some complications, the central (univariate) arguments in the proof of 
Theorem 14.21 can be directly generalised to multivariate polynomial rings. How- 
ever, this requires again certain assumptions on the commutation relations ^ in 
order to ensure that all needed computations are possible. 

Definition 4.3 The polynomial algebra of solvable type (P, -<) has centred com- 
mutation relations, if(i} there exists afield k C 7?. lying in the centre of TZ, (ii) the 
functions in ( 15a I ) are of the form p^{r) ~ p^{r)r with functions : 7^ — > k 
and ( Hi) we have r^^ G k in dJFJl. 

Using Konig's Lemma, Kredel proved in his thesis |42| the following version 
of Hilbert's Basis Theorem. 
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Theorem 4.4 Let [V, -<) be a polynomial algebra of solvable type with centred 
commutation relations over a left Noetherian coefficient ring TZ. Then V is left 
Noetherian, too. 

A third proof assumes that the ring V possesses a filtration S. Using an ap- 
proach detailed in |9 1 for the special case of the Weyl algebra (but which does not 
use any special properties of the Weyl algebra), one obtains the following general 
result where it is not even necessary to assume that is a polynomial ring. As 
discussed in Example 13.71 this result covers many of the polynomial algebras of 
solvable type that have appeared in the literature so far. 

Theorem 4.5 Let S be a filtration on the ring V. If the associated graded ring 
gTjjV is left Noetherian, then V is left Noetherian, too. 

Because of Condition (iii) in Definition 13. II we can define Grobner bases for 
ideals in algebras of solvable type. In the case that 7^ is a (commutative) field k, 
this is straightforward and from now on we will restrict to this case; the general 
case will be discussed only in SectionfTTl 

Definition 4.6 Let {V, ^) be a polynomial algebra of solvable type over afield k 
and I <Z V a left ideal. A finite set Q C V is a Grobner basis of X {for the term 
order <), if (le^ Q) — le^ X. 

For the ordinary multiplication this definition reduces to the classical one. The 
decisive point, explaining the conditions imposed in Definition l3.1l is that normal 
forms with respect to a finite set c "P may be computed in algebras of solvable 
type in precisely the same way as in the ordinary polynomial ring. Assume we 
are given a polynomial f E P such that le^ g \ le^ / for some g E G and set 
= Ic^ / — le^ g. If we consider — -k g, then by (iii) le^ — le^ /. 
Setting d = Ic^ //Ic^ 5^, we find by (ii) that le^ (/ — c?.g^) -< le^ /. Hence 
we may use the usual algorithms for computing normal form; in particular, they 
always terminate by the same argument as in the ordinary case. Note that in general 
d Ic^ // Ic^ 5, if r 7^ 1 in (|7]i, and that normal form computations are typically 
more expensive due to the appearance of the additional polynomial hin^. 

The classical Grobner basis theory can be extended straightforwardly to poly- 
nomial algebras of solvable type |f2l [T2l[T3ll4 1 Il42|[44, i45 1 . as most proofs are based 
on the computation of normal forms. The remaining arguments mostly take place 
in the monoid Nq and thus can be applied without changes. In particular, a triv- 
ial adaption of the standard approach leads to the following result crucial for the 
termination of Buchberger's algorithm. 

Theorem 4.7 Let {V, <) be a polynomial algebra of solvable type over afield. 
Then V is a left Noetherian ring and every left ideal X QV possesses a Grobner 
basis with respect to -<. 

Remark 4.8 Even in the case of a coefficient field we cannot generally expect V 
to be a right Noetherian ring, too; a concrete counterexample is provided again 
by McConnell and Robson pTl Example 1. 2. II]. In the proof of Theorem 14.71 
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one essentially uses that in normal form computations one always multiplies with 
elements of V from the left. Because of the already above mentioned left-right 
asymmetry of Definition 13.11 right ideals show in general a completely different 
behaviour. In order to obtain right Noetherian rings we must either adapt corre- 
spondingly our definition of a solvable algebra or impose additional conditions on 
the commutation relations (|5]l. 

The simplest possibility is to require that all the maps in (|5]l are automor- 
phisms (by Proposition l3.3l it suffices, if the maps pi in dSal l satisfy this condition). 
In this case we have k ★ + k = ★ k + k for all variables xi implying that we 
can rewrite any polynomial / c^x^^ in the "reverse" form / = ★ c^j. 

Now a straightforward adaption of the classical proof of Theorem 14.7 1 shows that 
the ring V is also right Noetherian. <\ 

We do not give more details on Grobner bases, as they can be found in the 
above cited references. Instead we will present in the next section a completely 
different approach leading to involutive bases. 



5 Involutive Bases 

We proceed to define involutive bases for left ideals in polynomial algebras of 
solvable type. In principle, we could at once consider submodules of free modules 
over such an algebra. As this only complicates the notation, we restrict to the ideal 
case and the extension to submodules goes as for Grobner bases. 

Definition 5.1 Let {V, -<) be a polynomial algebra of solvable type over a field k 
and 2 QV a non-zero left ideal. A finite subset Ti. C T is a weak involutive basis 
of X for an involutive division L on Ng , ;/ its leading exponents le^ H. form a 
weak involutive basis of the monoid ideal le^ X. The subset Ti is a (strong) invo- 
lutive basis of X, if le^ 7i is a strong involutive basis of le^ X and no two distinct 
elements of 7i have the same leading exponents. 

Remark 5.2 This definition of an involutive basis is different from the original one 
given by Gerdt and Blinkov [20 1 . Firstly, the notion of a weak basis is new. Sec- 
ondly, we do not require that an involutive basis is completely autoreduced. Finally, 
our approach is a natural extension of Definition l4.6l of a Grobner basis in the ring 
V whereas the approach of Gerdt and Blinkov [20] is closer to the constructive 
characterisation of Grobner bases via 5-polynomials. However, we will see below 
that essentially both approaches are equivalent. <\ 

Definition |5 . 1 [ implies immediately that any weak involutive basis is a Grobner 
basis. As in Section|2] we call any finite stt T <Z V (weakly) involutive, if it is a 
(weak) involutive basis of the ideal (JT) generated by it. 

Definition 5.3 Let J- (ZVX {0} be a finite set and L an involutive division on Nq . 
We assign to each element f <^ T a set of multiplicative variables 
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The involutive span of T is then the set 

(^)l,^ = E k[^L,^,^(/)] * / C (^) . (18) 

An important aspect of Grobner bases is the existence of standard represen- 
tations for ideal elements. For (weak) involutive bases a similar characterisation 
exists and in the case of strong bases we even obtain unique representations. 

Theorem 5.4 Let X Q V be a non-zero ideal, Ti <Z T\ {0} a finite set and L an 
involutive division on Nq. Then the following two statements are equivalent. 

(i) The set TL is a weak involutive basis of I with respect to L and -<. 

(ii) Every polynomial / G I can be written in the form 

f=Y,Ph^h (19) 

hen 

where the coefficients Ph G k[XL.-H,^(/i)] satisfy le^ {Ph * h) ^ Ic^ f for all 
polynomials h d Ti.. 

Ti. is a strong involutive basis, if and only if the representation H19\l is unique. 

Proof Let us first assume that the set H is a weak involutive basis. Take an arbi- 
trary polynomial / G T. According to Definition 15.11 its leading exponent Ic^ / 
lies in the involutive cone C^jc^ of at least one element h £ H. Let = 
Ic^ / — Ic^ h and set /i = / — cx^ -k h where the coefficient c G k is chosen such 
that the leading terms cancel. Obviously, /i G I and le^ /i -< le^ /. Iteration 
yields a sequence of polynomials fi G 2. After a finite number of steps we must 
reach Jn = 0, as the leading exponents are always decreasing and by assumption 
the leading exponent of any polynomial in X possesses an involutive divisor in 
Ic^ H. But this implies the existence of a representation of the form (fT9] l. 

Now assume that H is even a strong involutive basis and take an involutive 
standard representation ( fT9] l. By definition of a strong basis, there exists one and 
only one generator h G H such that le^ {Ph -k h) ^ le^ /. This fact determines 
uniquely le^ Ph- Applying the same argument to f — {It^ Ph)-*rh shows by recur- 
sion that the representation ( fT9l ) is indeed unique. 

For the converse note that (ii) trivially implies that Ic^ / G (Ic^ for 
any polynomial / G I. Thus Ic^ X C (Ic^ Ti)^,^. As the converse inclusion is 
obvious, we have in fact an equality and W is a weak involutive basis. 

Now let us assume that the set H is only a weak but not a strong involutive 
basis of X. This implies the existence of two generators /ii,/i2 G H such that 
C_L,ie^ •^(le-; C Cije^ ^(16^ hi). Hence we have Im^ /12 = Im^ (ex'' * hi) 
for suitably chosen c G kand/x G Nq. Consider the polynomial /i2 — cx''*/ii G X. 
If it vanishes, we have found a non-trivial involutive standard representation of 0. 
Otherwise an involutive standard representation /12 — cx^ -k hi = J^hen Ph'*^ h 
with Ph G ^[XLM,<{h)] exists. Setting P'f^ — Ph for all generators h ^ /ii,/i2 
and P,j^ = P/ij +cx^, P'f^^ ~ P/i^ — 1 yields again a non-trivial involutive standard 
representation = J^hei-c Pfi * ^- The existence of such a non-trivial representa- 
tion of immediately implies that iT% cannot be unique. Thus only for a strong 
involutive basis the involutive standard representation is always unique. □ 
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Corollary 5.5 Let the set Ti. be a weak involutive basis of the left ideal X Q V. 
Then {H)l,-< =1- 

Example 5.6 It is not true that any set T with {!F)l.^ = X is a weak involutive 
basis of the ideal I. Consider in the ordinary polynomial ring 'k[x, y] the ideal I 
generated by the two polynomials /i = and /2 = y"^ + . If we order the 
variables as xi ~ x and X2 = y, then the set T = {/i,/2} trivially satisfies 
(JT) = T, as with respect to the Janet division all variables are multiplicative 
for each generator. However, le^ J- — {[0,2]} does not generate le^ T, as obvi- 
ously [2, 0] G le^ T \ ({[0, 2]}). Thus !F is not a weak Janet basis (neither is the 
autoreduced set T' = {y^,x^}, as x^y ^ (jr')j^). <\ 

Proposition 5.7 Let X Q V be an ideal and Ti. C V a weak involutive basis of it 
for the involutive division L. Then there exists a subset H' ^ Ti. which is a strong 
involutive basis of X. 

Proof If the set le^ 7i is already a strong involutive basis of Ic^ X, we are done. 
Otherwise Ji contains polynomials h\, such that Ic^ h\ |l,1c^ u le^ Con- 
sider the subset H! = H\ {/i2}- As in the proof of Proposition l2 . 8 l one easily shows 
that le^ H' ~ le^ H \ {le^ /12} is still a weak involutive basis of le^ X and thus 
H' is still a weak involutive basis of X. After a finite number of such ehminations 
we must reach a strong involutive basis. □ 

Given this result, one may wonder why we have introduced the notion of a 
weak basis. The reason is that in more general situations like computations in 
local rings or polynomial algebras over coefficient rings (treated in later sections) 
strong bases rarely exist. 

Definition 5.8 Let J- a V be a finite set and L an involutive division. A polyno- 
mial g is involutively reducible with respect to T, if it contains a term x^ such 
that le^ / |L,ie^ r li-for some f ^ T. It is in involutive normal form with respect 
to J-, if it is not involutively reducible. The set T is involutively autoreduced, if no 
polynomial j T contains a term x^ such that another polynomial f E J-\ {/} 
exists with le^ /' |l,1o_; fJ- 

Remark 5.9 The definition of an involutively autoreduced set cannot be formulated 
more concisely by saying that each / e is in involutive normal form with respect 
to \ {/}. If we are not dealing with a global division, the removal of / from 
will generally change the assignment of the multiplicative indices and thus affect 
the involutive divisibility. <1 

An obstruction to involution is a polynomial g e (JT) \ possessing 
a (necessarily non-involutive) standard representation with respect to J^. We will 
later see that these elements make the difference between an involutive and an 
arbitrary Grobner basis. 

Example 5.10 Consider the set T ^ {/i,/2,/3} C ]k;[a::,y,2:] with the polynomials 
/i = — xy, f2 = yz — x and /a — y^ ~ z. For any degree compatible term order, 
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the leading terms of /2 and /a are unique. For /i we have two possibihties: if we 
use the degree lexicographic order (i. e. for x ^ y ^ z), it is z"^, for the degree 
inverse lexicographic order (i. e. for x y y y z) the leading term is xy. 

In the first case, (■^)./.^(jegiex ^ ^'^^ ^^^^ term order is a Janet 

basis, i. e. an involutive basis with respect to the Janet division, although we have 
not yet the necessary tools to prove this fact. In the second case, f4^ — z^~x'^ — 
zfi + xf2 G {J-) does not possess a standard representation and T is not even 
a Grobner basis. Adding to T yields a Grobner basis Q of {T), as one may 
easily check. But this makes z non-multiplicative for /2 and /s — z/2 is now an 
obstruction to involution of Q, as it is not involutively reducible with respect to the 
Janet division. In fact, the set ~ {/i, /2, /s, /4, /s} is the smallest Janet basis 
of I for this term order, as it is not possible to remove an element. Note that this 
second basis is not only larger but also contains polynomials of higher degree. <\ 

Remark 5.11 If is a Grobner basis of the ideal 2, then any element of I has a 
standard representation. But this does not imply that for a given division L the 
ideal I is free of obstructions to involution. In order to obtain at least a weak 
involutive basis, we must add further elements of X to ^ until (le_< G)l = le^ I- 
Obviously, this observation allows us to reduce the construction of a polynomial 
involutive basis to a Grobner basis computation plus a monomial completion. But 
we will see later that better possibilities exist. 

It follows that in general involutive bases are not reduced Grobner bases, as 
we already observed in Example 15. 101 For ^degiex the set was simultaneously 
a Janet basis and a reduced Grobner basis. But for ^deginviex the reduced Grobner 
basis is T U {/4}, whereas a Janet basis requires in addition the polynomial /s. 
We will see in Part II that this "redundancy" in involutive bases is the key for their 
use in the structure analysis of polynomial ideals and modules. <l 

It often suffices, if one does not consider all terms in g but only the leading 
term lt_< g: the polynomial g is involutively head reducible, if le^ / |l,ic^ r le^ g 
for some / G JT. Similarly, the set is involutively head autoreduced, if no leading 
exponent of an element / S is involutively divisible by the leading exponent of 
another element /' G \ {/}. Note that the definition of a strong involutive basis 
immediately implies that it is involutively head autoreduced. 

As involutive reducibility is a restriction of ordinary reducibility, involutive 
normal forms can be determined with trivial adaptions of the familiar algorithms. 
The termination follows by the same argument as usual, namely that any term 
order is a well-order. If g' is an involutive normal form of g ^ V with respect to 
the set for the division L, then we write g' = NFjr j^^^i^g), although involutive 
normal forms are in general not unique (like ordinary normal forms). Depending 
on the order in which reductions are applied different results are obtained. 

The ordinary normal form is unique, if and only if it is computed with re- 
spect to a Grobner basis; this property is often used as an alternative definition of 
Grobner bases. The situation is somewhat different for the involutive normal form. 

Lemma 5.12 The sum in ( l7iSI ) is direct, if and only if the finite set J-' C V \ {0} is 
involutively head autoreduced with respect to the involutive division L. 
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Proof One direction is obvious. For the converse, let /i, /2 be two distinct ele- 
ments of and Xi = their respective sets of multiplicative variables 
for the division L. Assume that two polynomials Pi £ k[Xi] exist with Pi * fi = 
P2 * /2 and hence le^ (Pi * /i) = Ic^ (P2 * /2)- As the multiplication ★ respects 
the term order -<, this implies that Cl,\c^ :f{\c^ fi) D Cl,\c^ ^^(10^ /2) ^ 0- Thus 
one of the involutive cones is completely contained in the other one and either 
le-; /i |l,1o^ le-; /2 or le^ /2 |l jc^ le^ A contradicting that T is involutively 
head autoreduced. □ 

Proposition 5.13 If the finite set T G V \ {0} is involutively head autoreduced, 
every polynomial g £ V has a unique involutive normal form NFjr.L.^(5)- 

Proof If is an involutive normal form of g, then obviously g G Con- 
versely, assume that g G {J-')l,^, i-c. the polynomial g can be written in the 
form g = X^/ej^^/ * / '^i'^^ Pf G ]k;[XijF,^(/)]. As is involutively head 
autoreduced, the leading terms of the summands never cancel (see the proof of 
Lemma l5.12b . Thus le^ g = le^ {Pf * /) for some f E T and any polynomial 
g G L.-< is involutively head reducible with respect to T. Each reduction step 
in an involutive normal form algorithm leads to a new polynomial g' G 
with le^ g' ^ le^ g. If the leading term is reduced, we even get le^ g' -< le^ g. As 
each terminating normal form algorithm must sooner or later reduce the leading 
term, we eventually obtain as unique involutive normal form of any g G {J-) 

Let gi and 52 be two involutive normal forms of the polynomial g. Obviously, 
51 — 52 G {^)l,-<- By definition of a normal form, neither gi nor 32 contain any 
term involutively reducible with respect to and the same holds for ,gi — (72 ■ Hence 
the difference gi — 92 is also in involutive normal form and by our considerations 
above we must have gi — 32 = 0. □ 

Proposition 5.14 The ordinary and the involutive normal form of any polynomial 
g £ V with respect to a finite weakly involutive set T <ZV \ {0} are identical. 

Proof Recalling the proof of the previous proposition, we see that we used the 
assumption that IF was involutively head autoreduced only for proving the exis- 
tence of a generator f E T such that le^ / |l,io^ r le^ g for every polynomial 
g G {J-)l,^- But obviously this property is also implied by the definition of a 
weak involutive basis. Thus by the same argument as above, we conclude that 
the involutive normal form with respect to a weakly involutive set is unique. For 
Grobner bases the uniqueness of the ordinary normal form is a classical property 
and any weak involutive basis is also a Grobner basis. As a polynomial in ordinary 
normal form with respect to F is trivially in involutive normal form with respect 
to JT, too, the two normal forms must coincide. □ 

Finally, we extend the notion of a minimal involutive basis from Nq to V. This 
is done in the same manner as in the theory of Grobner bases. 

Definition 5.15 Let T <Z V be a non-zero ideal and L an involutive division. An 
involutive basis Ti of T with respect to L is minimal, ;/ le^ H is the minimal 
involutive basis of the monoid ideal le_< T for the division L. 
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By Proposition 12. Ill we find that for a globally defined division like the Pom- 
maret division any involutive basis is minimal. Uniqueness requires two additional 
assumptions. First of all, our definition of an involutive basis requires only that it is 
involutively head autoreduced; for uniqueness we obviously need a full involutive 
autoreduction. Secondly, we must normalise the leading coefficients to one, i. e. 
we must take a monic basis. 

Proposition 5.16 Let I <Z V be a non-zero ideal and L an involutive division. 
Then I possesses at most one monic, involutively autoreduced, minimal involutive 
basis for the division L. 

Proof Assume that Tii and TL2 are two different monic, involutively autoreduced, 
minimal involutive bases of I with respect to L and -<. By definition of a min- 
imal involutive bases, this implies that Ic^TLi — As Hi and Ti,2 are 
not identical, we must have two polynomials hi G Hi and /i2 G H2 such that 
le^ hi = le^ h2 but hi ^ /i2. Now consider the polynomial h — hi — h2 G I. 
Its leading exponent must lie in the involutive span of le^ Hi — le^ H2- On the 
other hand, the term It^h must be contained in either hi 01 /i2. But this implies 
that either Hi or H2 is not involutively autoreduced. □ 

6 Monomial Completion 

We turn to the question of the actual construction of involutive bases. Unfortu- 
nately, for arbitrary involutive division no satisfying solution is known so far. In 
the monomial case, one may follow a brute force approach, namely performing 
a breadth first search through the tree of all possible completions. Obviously, it 
terminates only, if a finite basis exists. But for divisions satisfying some additional 
properties one can design a fairly efficient completion algorithm. 

The first problem in constructing an involutive completion of a finite subset 
Af C Wq for a division L is to check effectively whether J\f is already involutive. 
The trouble is that we do not know a priori where obstructions to involution might 
lie. If we denote by Ij G INq the multi index where all entries are zero except the 
jth one which is one, then the multi indices v + Ij with v G M and j G Nl,n'{v) 
are a natural first guess. 

Definition 6.1 The finite set M C Mg is locally involutive /or the involutive divi- 
sion L, if v +\j G {M) L for every non-multiplicative index j G Nl,N'{v) of every 
multi index v G M . 

Obviously, local involution is easy to check effectively. However, while (weak) 
involution obviously implies local involution, the converse does not necessarily 
hold. A concrete counter example was given by Gerdt and Blinkov ||20| . But they 
also discovered that for many divisions the converse is in fact true and thus for 
such divisions we can effectively decide involution. 

Definition 6.2 Let L be an involutive division and J\f C Ng a finite set. Let fur- 
thermore . . . , i^^'-*) be a finite sequence of elements of Af where every multi 
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index v^^^ with k < t has a non-multiplicative index jk G NL,jsf{v'^^^) such that 
v^^^ + Ij^. The division L is continuous, if any such sequence consists 
only of distinct elements, i. e. if v'^^^ ^ i^*-^-* for all k ^ L 

Proposition 6.3 For a continuous division L, any locally involutive set M C Ng 
is weakly involutive. 

Proof Let the set S contain those obstructions to involution that are of minimal 
length0 We claim that for a continuous division L all multi indices cr G are of 
the form v + Ij with v ^ M and j e ^lM^'A- This immediately implies our 
proposition: since for a locally involutive set all such multi indices are contained 
in {M) L, we must have = and thus (TV) = (A/')l- 

In order to prove our claim, we choose a <t E S for which no v E JV exists 
with a = ly + Ij. We collect in JVa all divisors ly e AA of cr of maximal length. 
Let be an element of Afc,; by assumption the multi index fi'^^^ — a ~ t/^^^ 
satisfies l/x^^^j > 1 and at least one non-multiplicative index ji G Nl,j^{i'^-^'>) 
exists with /ij^-* > 0. By the definition of S we have u^^^ + 1^^ G {■N')l- Thus a 
multi index i/^^' G M exists with v'^'^^ \L,j\f J^^^-* + Iji- This implies v'^'^^ \ a and 
we set /i(^) = (7 — v'^'^K By the definition of the set A/'ct we have < 
Hence i/^^) + l^- e {J\f) ^ for all j. 

Choose a non-multiplicative index ^2 G TVl jv'(i^'^') with /i^^-* > 0. Such an 
index exists as otherwise <j G {J^)l- By the same arguments as above, a multi 
index i/^^^ G TV exists with i/^^^ i'^'^^ + Ija and \ < |. We can iterate 
this process and produce an infinite sequence {v^^^ v'^'^\ . . . ) where each multi 
index satisfies u^^ G TV and |l,a/'J^^'^ + Ij. with j, G TVl^aa(i^^'^)- As TV 

is a finite set, the elements of the sequence cannot be all different. This contradicts 
our assumption that L is a continuous division: by taking a sufficiently large part 
of this sequence we obtain a finite sequence with all properties mentioned in Def- 
inition |62] but containing some identical elements. Hence a multi index G TV 
must exist such that a = v + \j. □ 

Lemma 6.4 The Janet and the Pommaret division are continuous. 

Proof Let TV C Nq be a finite set and (j^^*-*, . . . , i^^*-*) a finite sequence where 
+ with j G Nlm{^'^^ ) for 1 < i < t. 
We claim that for L = J, the Janet division, ^i^^ z^f*) implying that the 

sequence cannot contain any identical entries. Set k — max{i | /i^ Vi]. Then 
j < k, as otherwise j G TV,/^jv(i^(*+^^ ) entails j G Nj^j^{v'^'^^) contradicting our 
assumption that j is non-multiplicative for the multi index But j < kis also 
not possible, as then i^^'^^' < and so k cannot be multiplicative for 
There remains as only possibility j — k. In this case ^'j*^^'' — v^j^ + 1, as otherwise 
j could not be multiplicative for v^'^^^\ Thus we conclude that )^iex v''^^ 

and the Janet division is continuous. 

^ The length of a multi index v G Nq is the sum of its entries, i. e. the degree of the 
monomial x" . 
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The proof for the case L ~ P, the Pommaret division, is sHghtly more sub- 
tlel The condition j e iVp(i^(*)) implies that cls(i^(') + \j) = clsj/'^*' and if 
|p + Ij, then els > els i. e. the class of the elements of the 

sequence is monotonously increasing. If els = els i^^'^ = k, then the in- 

volutive divisibility requires that v'^^'^^ < v\f^ , i.e. among the elements of the 
sequence of the same class the corresponding entry is monotonously decreasing. 
And if finally v^^^^^ = v^f^K then we must have = + Ij, i. e. the length 

of the elements is strictly increasing. Hence all elements of the sequence are dif- 
ferent and the Pommaret division is continuous. □ 

Remark 6.5 In Remark IZ91 we discussed that for a global division a weak involu- 
tive basis of the sum Zi + T2 of two monoid ideals is obtained by simply taking 
the union of (weak) involutive bases of Ii and l2- As a more theoretical applica- 
tion of the concept of continuity, we prove now a similar statement for the product 
Ti ■ T2 and the intersection Ti n T2 in the special case of the Pommaret division. 
Let A/i be a (weak) Pommaret basis of Zi and J\f2 of T2. We claim that the set 
J\f = {fi + \ fi E Ml , V G 7V2} is a weak Pommaret basis of Ii ■ I2 and that the 
set M = {1cm (/i, J/) I /I G A/i , v G A/2} is a weak Pommaret basis of Ii 0X2- 

By Proposition l6.3l it suffices to show that the sets M and Sf, respectively, are 
locally involutive for the Pommaret division. Thus we take a generator /i + G N, 
where we assume for definiteness that els /i < els v, and a non-multiplicative index 
ji > els [fx^ v) — els /i of it. Then ji is also non-multiplicative for /i G A/i alone 
and the Pommaret basis A/i must contain a multi index jjS^^ which involutively 
divides + 1^^ . If we are lucky, then the generator /i*^^) + G AA is an involutive 
divisor of /i + + 1^^ , too, and we are done. 

Otherwise, there exists an index ki > els ly such that (/i — i-S^'')ki > 0. In this 
case the Pommaret basis A2 must contain a multi index z^'^^ which involutively 
divides i^+l^^. Again, if we are lucky, then /i'^^^+i^'^' G A/" is an involutive divisor 
of 11 + 1'+ Ij-^ and we are done. Otherwise, there are two possibilities. There could 
be an index j2 > els /x^^^ such that {fj, + i/+ Ij^ — /x'^' + i'''^^)j2 > entailing the 
existence of a further generator /i'^) G A/i which involutively divides /I'-^^ + lj^. Or 
there could exist an index ^2 > els v'^^^ such that {ji + v+lj-^ — fi^^'' +i>'^^^)k2 > 
implying that there is a multi index z^^^^ G A/2 involutively dividing z^*^^^ + l^^. 

Continuing in this manner, one easily sees that we build up two sequences 
(fi, . . . ) C A/i and (v, i/^^), i/^^), . . . ) C 7V2 as in the definition of a 

continuous division. Since both Pommaret bases are finite by definition and the 
Pommaret division is continuous by Lemma 16741 no sequence may become infinite 
and the above described process must stop with an involutive divisor of fi + iy+lj-^. 
Hence TV is locally involutive and a weak Pommaret basis of Zi • l2- The proof 



^ It is tempting to tackle the Pommaret division in the same manner as the Janet division 
using -<,evlex instead of -<iex; in fact, such a "proof" can be found in the literature. Unfor- 
tunately, it is not correct: if = i/'-'^ + Ij, then -irevlex z^'"' although the latter 
multi index is a divisor of the former one (^levlex is not a term order!). Thus the sequences 
considered in the application of Definition 16.21 to the Pommaret division are in general not 
strictly ascending with respect to -irevlex- 
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for J\f goes completely analogously replacing at appropriate places the sum of two 
multi indices by their least common multiple. <\ 

Definition 6.6 Let L be a continuous involutive division and TV C Ng a finite 
set of multi indices. Choose a multi index v £ M and a non-multiplicative index 
j € NL.j\f{^) such that: 

(ii) if there exists /i G A/" and k G Nj^j^fdi) such that /i + 1/j | i/ + Ij but 
H+lk^v + Ij, then fi+lkE {JV)l- 

The division L is constructive, if for any such set Af and any such multi index 

V + \j no multi index p £ (TV) l with i/ + Ij G CL,_^fi]{p}{p) exists. 

In words, constructivity may roughly be explained as follows. The conditions 
imposed on ly and j ensure a kind of minimality: no proper divisor of ly + Ij is 
of the form /i + 1^ for a /i G TV and not contained in the involutive span (TV)l. 
The conclusion implies that it is useless to add multi indices to TV that lie in some 
involutive cone, as none of them can be an involutive divisor of ly + Ij . An effi- 
cient completion algorithm for a constructive division should consider only non- 
multiplicative indices. 

Lemma 6.7 Any globally defined division ( and thus the Pommaret division ) is con- 
structive. The Janet division is constructive, too. 

Proof For a globally defined division the proof is very simple. For any multi index 
p G (TV) L there exists a multi index /i G TV such that p G Cl{ij-)- As for a globally 
defined division the multiplicative indices are independent of the reference set, 
we must have by the definition of an involutive division that Cl{p) Q Cl{ij)- 
Hence adding such a multi index to TV cannot change the involutive span and if 

V + \j ^ (TV) L, then also v + \j ^ (TV U {p})l- This implies constructivity. 

The proof of the constructivity of the Janet division is more involved. The basic 
idea is to show that if it was not constructive, it could not be continuous either. Let 
TV, V, j be as described in Definition 16.61 Assume for a contradiction that a multi 
index p G (TV)j exists with + Ij G C,7Xu{p}(p)- We write p — v'-^^ + p for 
a multi index i>'^^^ G TV with p G Cj^j\f{i'^^'>). As -I- Ij ^ {J^).j, we must have 
\p\ > 0. Set X = 1/ + Ij — p and let m, I be the maximal indices such that p„i > 
and A; > 0, respectively. 

We claim that j > maxj™, I}. Indeed, if j < m, then vii' < Vm and, by 
definition of the Janet division, this implies that m ^ Nj_j\f{v'^^^), a contradiction. 
Similarly, we cannot have j < I, as then / ^ -^,/.AAu{p}(p)- Finally, j — I is not 
possible. As we know already that j > m, we have in this case that pi = v'p = Vi 
for all i > 2 and pj < Vj. Hence j G TVj^y{pj(z^) and this implies furthermore 
J G -^,/.AAu{p}(p)' a contradiction. 

We construct a sequence as in Definition |62] of a continuous division. Choose 
an index ji with Aj^ > and ji G Nj^j^{v^^^). Such an index exists, as otherwise 
+ Ij e Cj,^f{v^^^^) C (TV)j. We write v + lj = {y^^^ + lj,)+p + \- Ij,. 
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Because of \fi\ > 0, the multi index v^^^ + Ij^ is a proper divisor of + Ij and 
according to our assumptions ly'^^^ G Af exists with t/'^^ + Ij^ G ^.^^(z^'^'). 

By the same arguments as above an index j2 G Nj^j^{v^'^^) must exist with 
(/i + A - > and a multi index j/^^) e TV with i/^^^ + Ij^ e C,7,^(j^^^^). 

Thus we can iterate and produce an infinite sequence . . . ) such that 

everywhere |j__v''^^''' + Iji with G Nj^j^{i>''^'>). By the continuity of the 

Janet division all members of the sequence must be different. However, every multi 
index v'^^'^ is a divisor of u + \j, so only finitely many of them can be different. 
Thus the sequence must terminate which only happens, if + \j G Cj^j^{v'^''''>) for 
some i contradicting our assumptions. □ 

We present now an algorithm for determining weak involutive completions 
of finite sets Af C Ng. As mentioned above, for arbitrary involutive divisions, 
nobody has so far been able to find a reasonable approach. But if we assume that 
the division is constructive, then a very simple completion algorithm exists, the 
basic ideas of which go back to Janet. 



Algorithm 2 Completion in (Nq , +) 

Require: a finite set H C Nq , an involutive division L 
Ensure: a weak involutive completion M of M 

\: N 

2: loop 

3: 5 ^ + 1, I y&M,j<=i NlM")^ I,M)l) 
4: if5 = 0tlien 
5: return A/^ 
6: else 

7: choose /j. £ S such that S does not contain a proper divisor of it 
8: Af^Afu{fi} 
9: end if 
10: end loop 



The strategy behind Algorithm |2] is fairly natural given the results above. It 
collects in a set S all obstructions to local involution. For a continuous division L, 
the set Af is weakly involutive, if and only if iS = 0. Furthermore, for a construc- 
tive division L it does not make sense to add elements of (AA) l to AA in order to 
complete. Thus we add in Line /8/ an element of S which is minimal in the sense 
that the set S does not contain a proper divisor of it. 

Proposition 6.8 Let the finite set Af C Ng possess a finite (weak) involutive com- 
pletion with respect to the constructive division L. Then Algorithm |2] terminates 
with a weak involutive completion Af of Af. 

Proof If Algorithm |2] terminates, its correctness is obvious under the made as- 
sumptions. The criterion for its termination, iS 0, is equivalent to local involu- 
tion of Af. By Proposition 16.31 local involution implies for a continuous division 
weak involution. Thus the result TV" is a weak involutive completion of TV, as by 
construction TV C Af C (Af). 
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If the input set Af is already involutive. Algorithm |2] leaves it unchanged and 
thus obviously terminates. Let us assume that TV is not yet involutive. In the first 
iteration of the loop a multi index of the form ji — v + Ij is added to N. It is 
not contained in {N)l and S does not contain a proper divisor of it. If Ml is an 
arbitrary involutive completion of TV, it must contain a multi index \ ^ J\f such 
that A a*- We claim that X = fi. 

Assume on the contrary that \ ^ [i. Since Ml C (TV), the multi index A 
must lie in the cone of a generator v^^^ G M . We will show that, because of the 
continuity of i, A G {M)l, contradicting the constructivity of L. If z^*^^^ \L,j\f\ 
we are done. Otherwise we write A = i^'^^ + p^^^ for some multi index p'^^^ G Nq . 
By construction, a non-multiplicative index ji G ^l,a/'(^^^'') with p)!'^ > must 
exist. Consider the multi index v^^^^ + Ij-^ . Because of z/*^^) + Ij^ | A, the multi index 
z/^^) + Ijj^ is a proper divisor of fj.. Since the set S does not contain any proper 
divisor of ^, we must have t^*^^^ + G (A/") l- Thus a multi index i/*^^^ GTVexists 
such that 

By iteration of this argument, we obtain a sequence v'^'^\ . . . ) where 

each element z^''^ G TV is a divisor of A and where zy('+i) v'^^^ + \j- with a 
non-multiplicative index ji G TVi_^(z/'^*'). This sequence cannot become infinite 
for a continuous division, as A possesses only finitely many different divisors and 
all the multi indices z^'*^ in arbitrary finite pieces of the sequence must be different. 
But the sequence will only stop, if some z/^'^ G TV exists such that z/^*' \l.i</ A and 
hence we must have that A G (TV) l ■ 

Thus every weak involutive completion TVl of the given set TV must contain 
the multi index v + Ij. In the next iteration of the loop. Algorithm [2] treats the 
enlarged set TVl = TV U {v + Ij}. It follows from our considerations above that 
any weak involutive completion Ml of TV is also a weak involutive completion of 
TVl and hence we may apply the same argument again. As a completion Ml is by 
definition a finite set, we must reach after a finite number k of iterations a weak 
involutive basis Mu of (TV). □ 

Note the crucial difference between this result and the termination proof of 
Buchberger's algorithm for the construction of Grobner bases. In the latter case, 
we can show the termination for arbitrary input, i. e. the theorem provides a con- 
structive proof for the existence of such a basis. Here we are only able to prove 
the termination under the assumption that a finite (weak) involutive basis exists; 
the existence has to be shown separately. For example. Lemma |2. 141 guarantees us 
that any monoid ideal possesses a finite weak Janet basis. 

Recall that by Proposition 12.81 anv weak involutive basis can be made strongly 
involutive by simply eliminating some redundant elements. Thus we obtain an 
algorithm for the construction of a strong involutive basis of (TV) by adding an in- 
volutive autoreduction as last step to Algorithm[2] Alternatively, we could perform 
the involutive autoreduction as first step. Indeed, if the input set TV is involutively 
autoreduced, then all intermediate sets TV constructed by Algorithm |2] are also in- 
volutively autoreduced. This is a simple consequence of the second condition in 
Definition 12. II of an involutive division that involutive cones may only shrink, if 
we add elements to the set TV. 
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Remark 6.9|1 While we just stated that it suffices to perform an involutive autore- 
duction as either first or last step in Algorithm |2l we now analyse for later use 
what happens, if we involutively autoreduce AT every time a new element has been 
added to it. The termination argument given in the proof of Proposition 16 . 8 1 does 
not remain valid after this modification and we must provide an alternative proof. 

Let again TVl = {m'^-* , ■ • ■ , Z^''^-' } be a weak involutive completion of the input 
set TV. If we denote by A/i the value of Af after the ith iteration of the loop, then 
it was shown in the proof of Proposition 16.81 that Ml is also a weak involutive 
completion of any set Afi. As by definition TVz, is finite and each Mi is a subset 
of it, the only possibility for non-termination is the appearance of a cycle, i. e. the 
existence of values fco. ^ such that Mk+e = Mk for all k > fco- 

Assume that in some iteration of the loop the multi index /i^''"' is added to 
M and that in the subsequent involutive autoreduction some elements of M are 
eliminated (in order to have a cycle this must indeed happen). The first step in the 
autoreduction must be that some multi index fi'-^^ is eliminated, because /i*^*^) is 
an involutive divisor of it. Indeed, by Condition (ii) in Definition 12.11 any other 
reduction would have been possible already before the insertion of /i'^*^-' and thus 
the previous involutive autoreduction would not have been finished. 

Since /i'^*^-' has been added to M, there must exist some multi index /i^^i) G M 
such that /i^*^^ = + p. Furthermore, we know that /i^^^ = n^'^'f + a for some 
multi index a with |a-| > and thus /i*^^^ = n^°-^'> + a with a = a + p and 
\a\ > 1. As we are in a cycle, the multi index p^^'' must have been added to M 
in a previous iteration of the loop, say when analysing Mi- Thus /i'^-' cannot be 
involutively divisible by p,^"-'-^ and we must have o-j^ > for a non-multiplicative 
index ji e Nj;^ j^.{p'^°-^^). It cannot be that p'^"'^^ + Ij^ = p'^^K as |(t| > 1, 
and therefore /i'^''^' + Ij^ is a proper divisor of p'-^^ Hence Mi must contain an 
involutive divisor p^"-^^ of p^"^^ + Ij-^, as otherwise this multi index would have 
been added to TV" instead of p^^\ 

Obviously, p^°-^'> \ p^'''> and, decomposing p^'''> — p^"-^") + n, we conclude 
by the same reasoning as above that tt^^ > for some non-multiplicative in- 
dex j2 e -^L.A/'i (m'"^'')- Iteration of this argument yields an infinite sequence 
. . .) as in Definition 16.21 of a continuous division. However, since 
i is a continuous division and Ml a finite set, we arrive at a contradiction. Thus 
even with involutive autoreductions after each step Algorithm|2]terminates. -O 

In some sense our description of Algorithm|2]is not complete, as we have not 
specified how one should choose the multi index p in Line /7/, if several choices 
are possible. One would expect that different involutive completions are obtained 
for different choices. However, an interesting aspect of our proof of Pi'oposition l6.8l 
is that it shows that this is not the case. The choice affects only the order in which 
multi indices are added but not which or how many multi indices are added during 
the completion. A simple method for choosing p consists of taking an arbitrary 
term order ~< (which also could be changed in each iteration of the loop) and 
setting p — min^ S. 



* The following considerations are joint work with Vladimir Gerdt. 
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Corollary 6.10 IfAlgorithm\2\terminates, its output M is independent of the man- 
ner in which fi is chosen. Furthermore, if A/l is any weak involutive completion 
of J\f with respect to the division L, then J\f C Ml- 

Proof Consider the set L(7V) of all weak involutive completions of TV with respect 
to the division L and define 

= fl Ml. (20) 

We claim that Algorithm |2] determines this set N independent of the used term 
order. Obviously, this implies our corollary. 

In the proof of Proposition 16 . 8 1 we showed that the multi indices added in Al- 
gorithm|2]are contained in every weak involutive completion of M. Thus all these 
multi indices are elements of N. As our algorithm terminates with a weak involu- 
tive completion, its output is also an element of L(A/') and hence must be N. □ 

Any monoid ideal in Nq has a unique minimal basis: take an arbitrary basis 
and eliminate all multi indices having a divisor in the basis. Obviously, these ehm- 
inations do not change the span and the result is a minimal basis. Similarly we 
have seen in Section|2]that if a monoid ideal I C JNg has a finite involutive basis 
for a given division L, then a unique minimal involutive basis exists. By the same 
argument as in the proof of Corollarv l6.10l it can easily be constructed by taking 
the unique minimal basis of I as input for Algorithmic 

7 Polynomial Completion 

An obvious way to compute an involutive basis for an ideal T in a polynomial 
algebra (T-",*, -<) of solvable type goes as follows: we determine first a Grobner 
basis QofI and then with Algorithmic an involutive completion of le^ Q. In fact, 
a similar method is proposed by Sturmfels and White |56| for the construction of 
Stanley decompositions (cf. Part II). However, we prefer to extend the ideas behind 
Algorithm|2]to a direct completion algorithm for polynomial ideals, as we believe 
that this approach is more efficient. 

First, we need two subalgorithms: involutive normal forms and involutive head 
autoreductions. The design of an algorithm NormalFormi.^((7, 7i) determining 
an involutive normal form of the polynomial g with respect to the finite setH C V 
is trivial. We may use the standard algorithm for normal forms in the Grobner 
theory, if we replace the ordinary divisibility by its involutive version. Obviously, 
this does not affect the termination. Actually, for our purposes it is not even nec- 
essary to compute a full normal form; we may stop as soon as we have obtained a 
polynomial that is not involutively head reducible. 

The design of an algorithm InvHeadAutoReduceL,^(.?^) for an involutive 
head autoreduction of a finite set T is also obvious. Again one may use the stan- 
dard algorithm with the ordinary divisibihty replaced by involutive divisibility. 



30 



Werner M. Seller 



Based on these two subalgorithms, we propose Algorithm [3] for the compu- 
tation of involutive bases in V. It follows the same strategy as the monomial al- 
gorithm. We multiply each generator by its non-multiplicative variables. Then we 
decide whether or not the result is already contained in the involutive span of the 
basis; if not, it is added. This decision is effectively made via an involutive nor- 
mal form computation: the involutive normal form of a polynomial is zero, if and 
only if the polynomial Ues in the involutive span. As our goal is a strong involutive 
basis, we take care that our set is always involutively head autoreduced. 



Algorithm 3 Completion in (P, -<) 

Require: a finite set JT c P, an involutive division L 

Ensure: an Involutive basis Ti. oil — (jT) with respect to L and -< 

I: H ^ InvHeadAutoReduceL,-< (.7^) 

2: loop 

3: 5 ^ {xj -k h \ h € Ti., Xj e XL,n,< (h), Xj-kh ^ {H)l,<) 
4: if5 = 0tiien 
5: return 7i 
6: else 

7: choose g & S such that le^ g = mln^ S 

8: 9 ^ NormalFormi,.-< ((7, H) 

9: H'^ InvHeadAutoReduceL,^(H U {g}) 

10: end if 

1 1 : end loop 



The manner in which we choose in Line 111 the next polynomial g to be treated 
(we briefly write min^ iS for the minimal leading exponent of an element of S) 
corresponds to the normal selection strategy in the theory of Grobner bases. There, 
this strategy is known to work very well for degree compatible term orders but not 
so well for other orders like the purely lexicographic one. Whereas for Grobner 
bases the selection strategy concerns only the efficiency of the computation, we 
will see below that here the use of this particular strategy is important for our ter- 
mination proof. With more refined and optimised versions of the basic completion 
Algorithm|3]one can circumvent this restriction lf5i lT4l[T9l . but we will not discuss 
this highly technical question here. 

Definition 7.1 A finite set T C P is locally involutive /or the division L, if for ev- 
ery polynomial f ^ T and for every non-multiplicative variable Xj £ 
the product Xj ★ / has an involutive standard representation with respect to J-. 

Note that for an involutively head autoreduced set JF, we may equivalently de- 
mand that Xj -k f £ because of Lemma 15. 12l this automatically implies 
the existence of an involutive standard representation. In fact, the criterion appears 
in this form in Line 131 of Algorithm |3] In any case, local involution may be ef- 
fectively verified by computing an involutive normal form of Xj * / in the usual 
manner, i. e. always performing head reductions. 
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Proposition 7.2 If the finite set J- <Z V is locally involutive for the continuous 
division L, then {J-)l.~^ = {^)- 

Proof We claim that if the set ^ is locally involutive (with respect to the con- 
tinuous division L), then every product x'* ★ /i of an arbitrary term a;^ with a 
polynomial fi E T possesses an involutive standard representation. This claim 
trivially entails our proposition, as any polynomial in [T) consists of a finite lin- 
ear combination of such products: adding the corresponding involutive standard 
representations shows that the polynomial is contained in (JT) 

In order to prove our claim, it suffices to show the existence of a representation 

* A = E * / + E ^-■/^^ * ^21) 

where Pf e k[XL,j^,M)] and le^ {Pf * f) = le^ {x'' * fi) (or Pf = 0) and 
where the coefficients c^j- G k vanish for all multi indices v G Nq such that 
le^ (x'' * f ) h Ic^ (a;^ * /i)- Our claim follows then by an obvious induction. 

If e ]k;[XLjF.^(/i)], i.e. it contains only variables that are multiplicative 
for le^ /i, nothing has to be shown. Otherwise we choose a non-multiplicative 
index ji G iVije_, jp-(le^ /i) such that /ij^ > 0. As is locally involutive, an 
involutive standard representation Xj-^ * fi = YlifeJ" -^/^' * ^ exists. Let T^Q T 
contain all polynomials /2 such that Ic^ (-f/a' * Ji) = Ic^ (xj^ ★ /i). If we have 
2;A'-iii g k[XL^;F_^(/2)] for all polynomials /2 G T^, then we are done, as at 
least Im^ (x^-i^i ^pf^) G k[Xi,^,^(/2)]. 

Otherwise we consider the subset T2 C T2 of polynomials for which 
x^^iji ^ k[XL,jc.^(/2)] and iterate over it. For each polynomial /2 G T'^ we 
choose a non-multiplicative index j2 G ^i.ic^ j^llc^ /2) such that (/^—Ij^ > 0. 
Again the local involution of the set T implies the existence of an involutive stan- 

dard representation Xj^ * f2 = 'l2fej=' P) * f- We collect in C all poly- 

(2) 

nomials such that le^ (-F/., * /s) — l*^^ i^j2 * /s)- If we introduce the multi 
index ly — le^ {xj^ * fi) — le^ f2, then le^ {x^^ -k fi) — le^ (x'^+^^^^i ~^J2 -k /s) 
for all /3 G T3. If x^+^-i^i -1^2 e k[XL,jF,^(/3)] for all /g G T3, we are done. 

Otherwise we continue in the same manner: we collect in a subset C 
all polynomials which are multiplied by non-multiplicative variables, for each 
of them we choose a non-multiplicative index G k[Xi,j(r,^(/3)] such that 
(/^ — Ijj — > 0, determine an involutive standard representation of xj^ ★ /s 

and analyse the leading terms. If they are still multiplied with non-multiplicative 
variables, this leads to sets JF4 C and so on. This process yields a whole tree 
of cases and each branch leads to a sequence (j/'-^-' — le^ /i, j/'-^-' = le^ /2, . . . ) 
where all contained multi indices j/^'^^ are elements of the finite set le^ J- and 
where to each j/*^''^ a non-multiplicative index jk G ^l,1c^ ) exists such that 

jyC^+i) 1^ jr jyC') + Ijj. . By the definition of a continuous division, this sequence 
cannot become infinite and thus each branch must terminate. But this implies that 
we may construct for each polynomial fi ^ T and each non-multiplicative vari- 
ables Xj G X[^^jr^^[fi) a representation of the claimed form (ISTT i. □ 
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Note that the proposition only asserts that the involutive span equals the normal 
span. It does not say that T is weakly involutive (indeed, the set T studied in Ex- 
ample |52] would be a simple counterexample). If 5 = X]/xgn" S/e.7=' ^t^-f^^ * / 
is an arbitrary polynomial in (JT), then adding the involutive standard represen- 
tations of all the products ★ / for which c^j ^ yields a representation 
g = J2feJ^-^f * / where Pf e iilXL.j^.^if)]- But in general it will not sat- 
isfy the condition le^ {Pf * /) ^ le^ g for all / e JT, as we cannot assume that 
we started with an ordinary standard representation of g. The satisfaction of this 
condition is guaranteed only for involutively head autoreduced sets, as there it is 
impossible that the leading terms cancel (Lemma |5.12| |. For such sets the above 
proof simplifies, as all the sets Ti consist of precisely one element and thus no 
branching is necessary. 

Corollary 7.3 For a continuous division L an involutively head autoreduced set 
J- C P is involutive, if and only if it is locally involutive. 

As in the proof of Proposition l6.8l local involution of H is obviously equivalent 
to the termination condition 5 = of the loop in Algorithm[3] Thus we are now 
in the position to prove the following result. 

Theorem 7.4 Let L be a constructive Noetherian involutive division and {P,-k, -<) 
a polynomial algebra of solvable type. Then Algorithm\3\tenninates for any finite 
input set J- with an involutive basis of the ideal X = {J-). 

Proof We begin by proving the correctness of the algorithm under the assumption 
that it terminates. The relation I — {Ti) remains valid throughout, although Ti, 
changes. But the only changes are the addition of further elements of 2 and invo- 
lutive head autoreductions; both operations do not affect the ideal generated by TL. 
When the algorithm terminates, we have 5 = and thus the output Ti is locally 
involutive and by Corollarv 17 . 31 involutive . 

There remains the problem of termination. Algorithm[3]produces a sequence 
(Til, ^^2, • • ■ ) with {Tii) — X. The set Tii+i is determined from Tii in Line 191. 
We distinguish two cases, namely whether or not during the computation of the 
involutive normal form in Line /8/ the leading exponent changes. If le^ g — le^ g, 
then (le^ 7ii) = (le^ as le^ g = le^ h + Ij for some h G Hi. Otherwise 

we claim that (le^ Hi) C (le^ Hi+i). 

By construction, g is in involutive normal form with respect to the set Hi im- 
plying that le^ g e (le^ Hi) \ (le^ Hi)L- If we had (le^ Hi) = (le^ Hi+i), a 
polynomial h G Hi would exist such that le^ g — le^ h + fi where the multi 
index /i has a non-vanishing entry fij for at least one non-multiplicative index 
j e Nl.\c^ Hi{h). This implies that le^ h + Ij ^ le^ g -< le^ g. But we choose 
the polynomial g in Line /7/ such that its leading exponent is minimal among all 
non-multiplicative products Xk*h with h E Hi; hence le^ g ^\e^h+lj. As this 
is a contradiction, we must have (le^ Hi) C (le^ Hi+i). 

So the loop of Algorithm [3] generates an ascending chain of monoid ideals 
(le^ Hi) C (le^ W2) C • • • C le^ X. As Nq is Noetherian, the chain must be- 
come stationary at some index N. It follows from the considerations above that in 
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all iterations of the loop after the iVth one le^ g = le^ g in Line /8/. At this stage 
Algorithm|3]reduces to an involutive completion of the monomial set le^ Hn us- 
ing Algorithmic — ^but with additional involutive autoreductions after each appear- 
ance of a new element. Indeed, in Line /7/ we choose the polynomial g such that 
le^ 5 is a possible choice for the multi index /i Algorithm|2]adds in Line /8/. Since 
we assume that our division is Noetherian, it follows now from Proposition 16.81 
together with Remark |6^ that Algorithmic] terminates (and our correctness proof 
above implies that in fact (le^ "Hn) = le^ T). □ 

Remark 7.5 If the division L is not Noetherian, then it may happen that, even when 
the ideal T = (JT) does possess a finite involutive basis with respect to L, Algo- 
rithm|3]does not terminate for the input We will see concrete examples for this 
phenomenon in Part II for the Pommaret division. 

The problem is that the existence of an involutive basis for le^ 2 does not 
imply that all subideals of it have also an involutive basis (as a trivial counter 
example consider (xy) C {xy, y^) with the Pommaret division). In such a case 
it may happen that at some stage of Algorithm |3] we encounter a basis Hi such 
that (le^ Hi) does not possess an involutive basis and then it is possible that the 
algorithm iterates endlessly in an attempt to complete Ic^ Hi- 

This observation entails that variations of Theorem l7.4l hold also for divisions 
which are not Noetherian. For example, we could assume instead that all subideals 
of le^ 2 possess an involutive basis. Alternatively, we could restrict to term orders 
of type uj. Then it suffices to assume that le^ 2 has an involutive basis. Indeed, 
now it is not possible that Algorithm|3]iterates endlessly within le^ Hi, as sooner 
or later an element g must be selected in Line /7/ with g ^ le^Hi- O 

Corollary 7.6 For a constructive Noetherian division L every ideal 2 V pos- 
sesses a finite involutive basis. 

Example 7. 7 Now we are finally in the position to prove the claims made in Ex- 
ample |5T0] With respect to the degree reverse lexicographic term order the Janet 
(and the Pommaret) division assigns the polynomial fi — — xy the multipUca- 
tive variables {x, y, z} and the polynomials /2 — yz — x and fs ^ y"^ — z the 
multiplicative variables {x, y}. Thus we must check the two non-multiplicative 
products: z/2 = yfi + xfs and z/3 = y/2 — /i. As both possess an involutive 
standard representation, the set S in Line /3/ of Algorithmic] is empty in the first 
iteration and thus is a Janet (and a Pommaret) basis of the ideal it generates. 

The situation changes, if we use the degree inverse lexicographic term or- 
der, as then It^ /i = xy. Now Xj^yr^^{fi) = {x}, Xj^jr^^{f.2) = {x,y,z} 
and X,/jF.^(/3) = {a;,?/}. In the first iteration we find S = {zfi}. Its involu- 
tive normal form is /a = z^ — x^ and we add this polynomial to T to obtain 
Hi = {/i, /2, /s, fi] (the involutive head autoreduction does not change the set). 
For /4 all variables are multiplicative; for the other generators there is one change: 
z is no longer multiplicative for /2. Thus in the second iteration S = {zf2}. It 
is easy to check that this polynomial is already in involutive normal form with 
respect to Hi and we obtain H2 by adding /s = yz"^ — xz to Hi. In the next iter- 
ation S is empty, so that H2 is indeed the Janet basis of (JF) for the degree inverse 
lexicographic term order. <1 
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Our proof of Theorem l7.4l has an interesting consequence which was first dis- 
covered by Apel f3l for the special case of the Pommaret division. Assume that 
the term order -< is of type u, i. e. for any two multi indices fi, v with < v only 
finitely many multi indices p''^ exist with /i -< p^^^ -< p^^^ -< • • • ^ v. Then 
even if our algorithm does not terminate, it determines in a finite number of steps 
a Grobner basis of the ideal T. 

Proposition 7.8 Let the term order -< be of type lo. Then Algorithm\3\determines 
for any finite input set T <Z V in a finite number of steps a Grobner basis of the 
ideal I = {J-). 

Proof Above we introduced the set Hn such that (le^ Ti.N+e) — (le^ "Hn) for all 
£ > 0. We claim that Hn is a Grobner basis of T. 

Let / G Z be an arbitrary element of the ideal. As 7i at is a basis of 2, we find 
for each h ^TLn polynomial g/, G V such that 

9h*h. (22) 

heT-lN 

Hn is a Grobner basis, if and only if we can choose the coefficients gh such that 
le^ {gh * /i) ^ Ic^ /. Assume that for / no such standard representation exists 
and let fj, = maxheHN {l^-; 9h + le^ >- le^ /. If we denote by Hn the set 
of all polynomials h ^ Hn for which le^g^ + le^ h = ji, then the identity 
TliheHN *h)=Q must hold and hence Hn contains at least two elements. 

For each element h G Hn we have /i G C(le^ h). As by construction the set Hn 
is involutively head autoreduced, the involutive cones of the leading exponents do 
not intersect and there must be at least one generator h G Hn such that some 
non-multiplicative variable Xj G Xl.-h„(/i) divides It^ g^. 

As -< is of type uj, after a finite number of steps the non-multiplicative product 
Xj -khis analysed in Algorithm[3] Thus for some ni > the set Hn+ui contains 
an element h' with le^ h! — le^ [xj * h). Let p, = le^ gj^, x^^^^ -kXj = cx^ + ri 
and h' = dxj -kh + r2 - Then we may rewrite 



^ I gj, 

9h 



cd 



•k [h' — r2) — dri -kh + [gj^ ~ Im^ gj^ ★ h . (23) 



As h' was determined via an involutive normal form computation applied to the 
product Xj -k h and as we know that at this stage of the algorithm the leading 
exponent does not change during the computation, the leading exponent on the 
right hand side of (l23l l is le^ {x^^^^ -kh'). If the term x^^^^ contains a non- 
multiplicative variable Xk G ^l,Wjv+„^ (^')' '^^ repeat the argument obtaining 
a polynomial h" G HN+ni+n2 such that le^ h" — le^ {xk * h'). 

Obviously, this process terminates after a finite number of steps, even if we do 
it for each h G Hn- Thus after £ further iterations we obtain a set Hn+i such 
that, after applying all the found relations (|23] |. / can be expressed in the form 
/ = J2heHN+e9h * h where still = max,,g-H„+j {le^ ff/t + le^ h}. Denote 
again by Hn+i ^ T^N+e the set of all polynomials h achieving this maximum. 
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By construction, no term It^ g,- with h e Hn+i contains a variable that is 
non-multiplicative for h. Thus we must now have /i G Cic^ (le^ h) for 

each h G 'Hn+i implying that Hn+i contains at most one element. But then it is 
not possible that fi >~ le^ /. Hence each polynomial f € V possesses a standard 
representation already with respect to Hn and this set is a Grobner basis. □ 

Note that in the given form this result is only of theoretical interest, as in gen- 
eral no efficient method exists for checking whether the current basis is already a 
Grobner basis. Using standard criteria would destroy all potential advantages of 
the involutive algorithm. For the special case of Pommaret bases, Apel 131 found a 
simple criterion that allows us to use a variant of Algorithm[3]for the construction 
of Grobner bases independent of the existence of a finite involutive basis. 

In contrast to the monomial case, one does not automatically obtain a minimal 
involutive basis by making some minor modifications of Algorithm[3] In particu- 
lar, it does not suffice to apply it to a minimal basis in the ordinary sense. Gerdt 
and Blinkov II2TI presented an algorithm that always returns a minimal involutive 
basis provided a finite involutive basis exists. While it still follows the same basic 
strategy of study all products with non-multiplicative variables, it requires a more 
complicated organisation of the algorithm. We omit here the details. 



8 Right and Two-Sided Bases 

We now briefly discuss the relation between left and right involutive bases and the 
computation of bases for two-sided ideals. We use in this section the following 
notations: the left ideal generated by c "P is denoted by (jF)'-'\ the right ideal 
by and the two-sided ideal by {(T)) and corresponding notations for the 

left, right and two-sided involutive span. 

Recall from Remark l4~8] that even with a coefficient field k it is not guaranteed 
that V is also right Noetherian and hence generally the existence of right Grobner 
bases for right ideals is not clear. However, we also noted that the ring V is always 
right Noetherian, if we assume that the maps pi : k — > k in dSal i are automor- 
phisms. In the sequel of this section we will always make this assumption. 

From a computational point of view, the theory of right ideals is almost iden- 
tical to the corresponding theory for left ideals. The left-right asymmetry in our 
definition of polynomial algebras of solvable type leads only to one complication. 
Suppose that we want to perform a right reduction of a term ax"^ with respect 
to another term cx^ with /i | 1/. This requires to find a coefficient 6 e k such 
that Ic^ (cx^ * bx'^~'^) — cpfj_{b)rf^^^-^ = a. Since, according to the above made 
assumption, all the maps are automorphisms, such a b always exists. 

It turns out |42 Sect. 4.11] that under the made assumption the results of 
Kandry-Rodi and Weispfenning [41, Sects. 4/5] remain valid for our larger class 
of non-commutative algebras and can be straightforwardly extended to involutive 
bases. For this reason, we will only discuss the case of involutive bases and do not 
treat separately Grobner bases. 
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Lemma 8.1 Let (P,*, ^) be an arbitrary polynomial algebra of solvable type 
where all the maps appearing in the commutation relations ( 15a I ) are automor- 
phisms. A polynomial f ^ V is (involutively) left reducible modulo a finite set 
T <Z V (with respect to an involutive division L), if and only if it is (involutively) 
right reducible (with respect to L). 

Proof Because of the made assumptions on the maps p^, reducibihty depends 
solely on the leading exponents. □ 

Proposition 8.2 Let Tii be a monic, involutively left autoreduced, minimal left in- 
volutive set and Tir a monic, involutively right autoreduced, minimal right involu- 
tive set for an involutive division L. If (Ti;)'^'-' = {TLr)^'^^ — T, then Tii = 7ir- 

Proof By definition of a minimal basis, the sets le^ Hi and le^ Hr are both min- 
imal involutive bases of the monoid ideal le^ 2 and thus are identical. Assume 
that (Til \ Hr) U {Hr \ Hi) ^ % and let / be an element of this set with minimal 
leading exponent with respect to -<. Without loss of generality, we assume that 
f eHi\ Hr- Because of the condition = (HrY''\ we have / e {Hr)'-[^^. 

Thus the (by Proposition 15 . 1 3 1 unique) right involutive normal form of / with re- 
spect to Hr is 0. This impUes in particular that / is right involutively reducible 
with respect to some h G Hr with le^ h ^ le^ /. 

If le^ h -< le^ /, then h £ Hi, too, as / was chosen as a minimal element of 
the symmetric difference of Hi and Hr- Hence, by Lemma ISTl / is also left invo- 
lutively reducible with respect to h (because of \e^Hi — le^ Hr the multiplicative 
variables of h are the same in both cases). But this contradicts the assumption that 
Hi is involutively left autoreduced. 

If le^ h = f = p, then we consider the difference g = f — h E T: both the 
left involutive normal form of g with respect to Hi and the right involutive normal 
form with respect to Hr must vanish. By construction, le^ g ^ 1^ and supp g C 
(supp/ U supp/i) \ {/i}. Since both Hi and Hr are assumed to be involutively 
autoreduced, no term in this set is involutively reducible by Ic^ Hi = Ic^ Hr and 
hence we must have supp g — 0, i. e. g = 0, a contradiction. □ 

A direct derivation of a theory of two-sided involutive bases along the lines of 
Section |5] fails, as two-sided linear combinations are rather unwieldy objects. A 
general polynomial / £ {{H)) for some finite setH C V is of the form 

"h 

/= ;^^^,*;i*r, (24) 

hen i=i 

with polynomials £ P, i. e. we must allow several summands with the same 
generator h. The definition of a unique involutive standard representation would 
require control over the numbers Uh which seems rather difficult. Therefore we 
will take another approach and construct left involutive bases even for two-sided 
ideals. The following proposition is an involutive version of Theorem 5.4 in [41]. 

Proposition 8.3 Let H C ("P,*, ~<) be a finite set and L an involutive division. 
Then the following five statements are equivalent. 
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(i) Ti. is a left involutive basis and (TY)^'^ = ((7Y)). 

(ii) TL is a right involutive basis and (H) = ((^))- 

(iii) Ti is a left involutive basis of {TL)^'-^ andbothh-kxi G {TLY'-^ andh-kc G (Ti)^'-* 
for all generators h E Ti., all variables Xi and all coefficients c S k. 

(iv) Ti, is a right involutive basis of {Ti) and both Xi-k h E (Ti) and c* h E 
(Ti) for all generators h E Ti, all variables Xi and all coefficients c S k. 

(v) A unique generator h E Ti. exists for every polynomial f E {{Ti.}) such that 
le^ h \l,ic-, n le^ /. 

Proof We begin with the equivalence of the first two statements, (i) implies that 
(n)'-l\ = (ny^^ = ({n)) and hence trivially C The same argu- 

ment as in the proof of Proposition l8.2l shows that we have in fact an equality and 
thus {'hCj^l^^ = {Ti)^^^ — {{H)), i. e. (ii). The converse goes analogously. 

Next we consider the equivalence of (i) and (iii); the equivalence of (ii) and 
(iv) follows by the same argument, (iii) is a trivial consequence of (i). For the 
converse, we note that (iii) implies that / ★ (c<) E {H)'-'-'' for all / E {'H)^^\ all 
terms t eT and all constants c E k. Indeed, as in the proof of Propo sition 1 3 . 3 1 we 
may rewrite the monomial ct as a polynomial in the "terms" x^^ ★ x*^ ★ • • ■ * a:*' 
with «i < 42 < ■ ■ ■ < and then apply repeatedly our assumptions. Obviously, 
this entails (i). 

The equivalence of (i) or (ii), respectively, with (v) is a trivial consequence of 
the definition of an involutive basis. □ 

We would Uke to exploit Statement (iii) for the construction of a left involu- 
tive basis for the two-sided ideal {{J-)). However, if the field k is infinite, then 
it contains an infinite number of conditions. In the sequel we will follow 1421 
Sect. 4.1 1] and make one further assumption about the polynomial algebra V. Let 
ko = {c G k I V/ G P : c ★ / = / ★ c} be the constant part of the centre of V. 

Lemma 8.4 ko is a subfield of k. 

Proof It is obvious that ko is a subring. Thus there only remains to show that with 
c e kp we have E ko, too. If c G ko, then Xi-k c = cxi, i. e. pi{c) = c and 
hi{c) = 0, for all 1 < i < n. Now on one hand Xi ★ (c~^ -k c) = Xi and on the 
other hand 

{xt ★ c"-^) ★ c = pi(c"^)pi(c)xi + chi{c^^) (25) 

{hi{c~^) kr c = chi{c~^) since c G ko). The associativity of ★ implies now that 
Pi{c~^) = and hi{c~^) — 0. Hence c^^ commutes with all variables Xi and 
it is easy to see that this entails c^^ G ko. □ 

We make now the assumption that either k^ = {ci, . . . , ci} is finite or that 
the extension k/ko is finite, i. e. that k is a finite-dimensional vector space over 
ko with basis {ci, . . . , c^}. In the latter case, it is easy to see that it suffices in (iii) 
to require that only all products h k cj lie in (H) as for c = ^j'^j ^i'^h 

Xj E ko we have h-k c ~ Yfj=i ^ji.^ * ^j)- 

These considerations lead to the simple Algorithm |4] below. It first constructs 
in Line /!/ a left involutive basis H. of the left ideal {T) (using Algorithm[3]). The 
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while loop in Lines /2-19/ extends the set to a left generating set of the two- 
sided ideal {{T)) according to our simplified version of statement (iii) in Propo- 
sition [O] Finally, we complete in Line /20/ this set to an involutive basis. Note 
that in Line /I/ it is not really necessary to compute a left involutive basis; any left 
Grobner basis would suffice as well. Similarly, an ordinary left normal form could 
be used in Lines /6/ and /12/, respectively; the use of InvLef tNormalFormi^^ 
anticipates the final involutive basis computation in Line /20/. 



Algorithm 4 Left Involutive basis for two-sided ideal in {V, *, ^) 
Require: finite set CV, involutive division L 
Ensure: left Involutive basis Ti. of {{J-}) 

1: ^ LeftInvBasis_L,^(J'^); S^H 

2: wliile 5 / do 



3: 

4: for all / G 5 do 

5: for i from 1 to n do 

6: h ^ InvLef tNormalFormi,^^ (/ * aji, 7i) 

7: if /i / then 

8: H^Hu{/i}; T^T\j{h} 

9: end if 

10: end for 

1 1 : for j from 1 to ^ do 

12: h <— InvLef tNormalForm_L,^ (/ * Cj, 7i) 

13: if /i 7^0 then 

14: H^nuih}; T ^ T VJ {h} 

15: end if 

16: end for 

17: end for 

18: S^T 



19: end while 

20: return LeftInvBasisL,^(7i) 



The termination of the while loop follows from the fact that under the made 
assumptions V is Noetherian and hence a finite generating set of {{J-Jj exists. In 
principle, we perform here a simple breadth-first search for it. The termination of 
the involutive bases computations in Lines III and /20/, respectively, depends on 
the conditions discussed in the last section. Thus the termination is guaranteed, if 
the division L is constructive and Noetherian. 



9 Involutive Bases for Semigroup Orders 

For a number of applications it is of interest to compute involutive or Grobner 
bases with respect to more general orders, namely semigroup orders (see Ap- 
pendix |A]l. This generalisation does not affect the basic properties of polynomial 
algebras of solvable type as discussed in Sect. [3] but if 1 is no longer the smallest 
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term, then normal form computations do no longer terminate for all inputs. So we 
can no longer apply Algorithm[3]directly for the determination of involutive bases. 

Example 9.1 The Weyl algebra W„ is the polynomial algebra in the 2n variables 
xi , . . . , and di, . . . ,dn with the following non-commutative product for aU 
1 < « < n we have di-kXi — Xidi + 1 and ★ is the normal commutative product in 
all other cases. It is easy to see that W„ is a polynomial algebra of solvable type 
for any monoid order For semigroup orders compatibility requires that 1 -< Xidi 
for all i. In |54| such orders are called multiplicative monomial orders. 

An important class of semigroup orders is defined via real weight vectors. Let 
C) G IR" IR" be such that G IR" is non-negative and let ^ be an arbitrary 
monoid order. Then we define x'^ 9" "^(C.C) if either/i-^ + i^-C < cr-^ + r-C 

or/i-^ + j/-C = o'-^ + r- C and x^d"^ -< x'^d'^ . This yields a monoid order, if and 
only if both ^ and C, are non-negative. A special case are the orders with weight 
vectors — ^) arising from the action of the algebraic torus (k*)" on the Weyl 
algebra. They have numerous applications in the theory of P-modules ||54]| . < 

As normal form computations do not necessarily terminate for semigroup or- 
ders, we must slightly modify our definitions of (weak) involutive or Grobner 
bases. The proof of Theorem 15 .41 (and consequently also the one of Corollarv l5.5l 
showing that a weak involutive basis of an ideal I is indeed a basis of I) requires 
normal form computations and thus this theorem is no longer valid. The same prob- 
lem occurs for Grobner bases. Therefore we must explicitly include this condition 
in our definition. 

Definition 9.2 Let (V, -k, ~<) be a polynomial algebra of solvable type where -< is 
an arbitrary semigroup order Let furthermore I Q V be a left ideal. A Grobner 
basis of T is a finite set Q such that (Q) — X and (le^ G) — le^ X. The set Q 
is a weak involutive basis of X for the involutive division L, if in addition the set 
le^ Q is weakly involutive for L. It is a (strong) involutive basis, if it is furthermore 
involutively head autoreduced. 

In the case of Grobner bases, a classical trick due to Lazard (43) consists of ho- 
mogenising the input and lifting the semigroup order to a monoid order on the ho- 
mogenised terms. One can show that computing first a Grobner basis for the ideal 
spanned by the homogenised input and then dehomogenising yields a Grobner ba- 
sis with respect to the semigroup order. Note, however, that in general we cannot 
expect that reduced Grobner bases exist. 

We extend now this approach to involutive bases. Here we encounter the ad- 
ditional difficulty that we must lift not only the order but also the used involutive 
division. In particular, we must show that properties like Noetherity or continuity 
are preserved by the lift which is non-trivial. For the special case of involutive 
bases in the Weyl algebra, this problem was first solved in |35|. 

Let (T',*, -<) be a polynomial algebra of solvable type where -< is any semi- 
group order that respects the multiplication We set V — h[xo, xi, . . . , Xn] and 
extend the multiplication Tkr to "P by defining that xq commutes with all other vari- 
ables and the elements of the field k. For a polynomial f ~ J2 C/^a^^ G "P of 
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degree q, we introduce as usual its homogenisation f^^^ = 

Conversely, for a polynomial / S "P we denote its projection to V as f = f\.xo=i- 

We denote by T the set of terms in V; obviously, it is as monoid isomorphic 
to N^+^ We use in the sequel the following convention. Multi indices in Nq 
always carry a tilde: /i = [/io, • • • , M"]- Ths projection to INq defined by dropping 
the first entry (i. e. the exponent of the homogenisation variable xo) is signalled by 
omitting the tilde; thus ^ = [/ii, . . . , For subsets J\f C Nq"*"^ we also simply 
write A/" = eJV} C N^J. 

We lift the semigroup order -< on T to a monoid order on T by defining 
xi^ -<h x^-, if either |/t| < |z>| or both \ fb\ = \i>\ and x^^ -< x^. It is trivial to check 
that this yields indeed a monoid order and that ("P, t^, -<h) is again a polynomial 
algebra of solvable type. For lifting the involutive division, we proceed somewhat 
similarly to the definition of the Janet division: the homogenisation variable xq is 
multipUcative only for terms which have maximal degree in xq. 

Proposition 9.3 Let L be an arbitrary involutive division on ^Nq . For any finite set 

M C Ng^^ and every multi index ji G J\f, we define Nj^ jv-(A) by- 

- e Nijij-{fi.), if and only if /xq = max-gj^{i^o}, 
-0 < i e if and only if i G Nl,m{ij)- 

This determines an involutive division L on Nq"*"^. 

Proof Both conditions for an involutive division are easily verified. For the first 
one, let p G Ci fl Ci /j-iiy') with e. J\f. If po = Mo = J^o. the first entry 
can be ignored, and the properties of the involutive division L implies the desired 
result. If po = Ho > vq, the index must be multiplicative for D contradicting 
po > fo- If Po is greater than both pQ and uq, the index must be multipUcative 
for both implying p^) = Uf). In this case we may again ignore the first entry and 
invoke the properties of L. 

For the second condition we note that whether a multiplicative index i > be- 
comes non-multiplicative for some element v £ M after adding a new multi index 
to N is independent of the first entry and thus only determined by the involutive 
division L. If the new multi index has a higher first entry than all elements of M, 
then becomes non-multipUcative for all elements in J\f but this is permitted. □ 

Now we check to what extent the properties of L are inherited by the lifted 
division L. Given the similarity of the definition of L and the Janet division, it is 
not surprising that we may reuse many ideas from proofs for the latter. 

Proposition 9.4 If L is a Noetherian division, then so is L. 

Proof Let A/" C Nq"'"^ be an arbitrary finite subset. In order to prove the existence 
of an Z-completion of A/', we first take a finite i-completion N c A'S of TV which 
always exists, as by assumption L is Noetherian. Next, we define a finite subset 
M' C {M) by setting 

M' = \p& N[J+^ I M e A/" A Mo < maxz/Q | n (TV) . 
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We claim that this set Af' is an i-completion of Af. By construction, we have both 
M' C (TV) and M C TV"', so that we must only show that M' is involutive. 

Let jl e (A/"') be arbitrary. By construction of TV"', we can find v G N' with 
v\j^ jjL. Moreover, the definition of Sf' guarantees that we can choose v in such 
a way that either vq = fiQ or = max-gj^, po < /io holds. In the former case, 
we trivially have i^li ji/, fj.', in the latter case we have G f^fiv) (see the proof 
of Proposition |93]l. Thus in either case jl G {Af') i- □ 

Proposition 9.5 If L is a continuous division, then so is L. 

Proof Let {i>^^^ , • • ■ , t'''^-') with t^'*) G TV" be a finite sequence as described in the 
definition of continuity. We first note that the integer sequence (j/q^"* , . . . , t'o'^'') is 
monotonically increasing. If z/q*^ is not maximal among the entries fiQ for jl G M, 
no multiplicative divisor of i/^'^ + Ij in J\f can have a smaller first entry: if i/q'^ 
is maximal, the index is multiplicative for and any involutive divisor in TV" 
must also be maximal in the zero entry. Thus it suffices to look at those parts of the 
sequence where equality in the zero entries holds. But there the inequality of the 
multi indices follows from the continuity of the underlying division L. □ 

Unfortunately, it is much harder to show that constructivity is preserved. We 
will do this only for globally defined divisions and the Janet division. 

Proposition 9.6 If the continuous division L is either globally defined or the Janet 
division, then the lifted division L is constructive. 

Proof We give a proof only for the case of a globally defined division. For the 
Janet division J one must only make a few modifications of the proof that J itself 
is constructive. We omit the details; they can be found in 1351 . 

We select a finite set A/" C Nq^^, a multi index // G TV" and a non-multiplicative 
index i of fi such that the conditions in the definition of constructivity are fulfilled. 
Assume that there exists a p G TV such that fi+li = p+a+r with p+a G ji/{p) 
and p+(T + f G a"u{p+ct} + ''')■ ^et L be a globally defined division. If i = 0, 
then /io + 1 = Po + (Tq + tq implies that (Tq = tq ^ 0: for (To > 0, we would 
have (0 is multiplicative for p) po > po ^ Po + <^o > Po- For do = and tq > 
a similar contradiction appears. If i > 0, the argumentation is simple. A global 
division is always constructive, as adding further elements to TV does not change 
the multiplicative indices. But the same holds for the indices fc > in the lifted 
division L. Thus under the above conditions p, + li G (TV)^ contradicting the 
made assumptions. □ 

Based on these results. Algorithm [3] can be extended to semigroup orders. 
Given a finite set g P, we first determine its homogenisation JTC*) e V and then 
compute an involutive basis of (i^^'*^) with respect to L and What remains to 
be done is first to show that the existence of a finite involutive basis is preserved 
under the lifting to P and then to study the properties of the dehomogenisation of 
this basis. 
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Proposition 9.7 If the left ideal I — {J-) C V possesses an involutive basis with 
respect to the Noetherian division L and the semigroup order ■<, then the left ideal 
X = ( J-"^''^ ) C ■p generated by the homogenisations of the elements in the finite 
set T possesses an involutive basis with respect to the lifted division L and the 
monoid order -<^. 

Proof By Theorem l4.7l the ideal I QV possesses a Grobner basis Q with respect 
to the monoid order ^/j. By Proposition 19.41 a finite L-completion M of the set 
le^;, Q exists. Moreover, as is a Grobner basis of I, the monoid ideals (le^^ Q) 
and le^^ I coincide. Thus M is an involutive basis of le_<^ X with respect to the 
lifted division L and an involutive basis H of T with respect to the division L is 
given by 

{x^'irglg eg A le^Ax~^*9) ^J^} ■ (26) 
This set is obviously finite. □ 

Hence our lifting leads to a situation where we can apply Theorem l7.4l Unfor- 
tunately, the dehomogenisation of the strong involutive basis computed in V does 
not necessarily lead to a strong involutive basis in V, but we obtain always at least 
a weak involutive basis and thus in particular a Grobner basis. 

Theorem 9.8 Let Ti. be a strong involutive basis of the left ideal Z C "P with 
respect to L and -<h- Then the dehomogenisation H. is a weak involutive basis of 
the left ideal I (-V with respect to L and -<. 

Proof For any / G Z an integer A; > exists such that / — Xq/^^'^ £ X. The 
polynomial / possesses a unique involutive standard representation 

~f=T.Phh (27) 

hen 

with e I'^t^L.ie^^ n(^y\ ™d le^^ {Pjji) <h le^^ /. Setting xq = 1 in (|27]i 
yields a representation of / with respect to the dehomogenised basi^ H. of the 
form / = X^/iew where Ph & ^[XL,ic^nih)] by the definition of the lifted 
division L. This obviously implies that (H) = I. By the definition of the lifted 
order and the homogeneity of the lifted polynomials, we have furthermore that 
le^ {Phh) :< le^ / and hence that Ic^ His a weak involutive basis of Ic^ X. Since 
all conditions of Definition |9.2| are satisfied, the set H is therefore indeed a weak 
involutive basis of the ideal X. □ 

Remark 9.9 For the Pommaret division P the situation is considerably simpler. 
There is no need to define a lifted division P according to Proposition l9.3l Instead 
we renumber xo to Xn+i and then use the standard Pommaret division on Nq^^. 
This approach implies that for all multi indices fl e Nq^^ with /i 7^ the equality 
Np (//) = Np (/i) holds, as obviously n + 1 is multiplicative only for multi indices 

' Note that the dehomogenised basis 7i is in general smaller than 7i, as some elements 
of may differ only in powers of xo- 
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of the form p, = £n+i, i- e. for which /i = 0. One easily sees that the above proof 
of Theorem l9.8l is not affected by this change of the division used in Nq^^ and 
hence remains true. < 

It is not a shortcoming of our proof that in general we do not get a strong 
involutive basis, but actually some ideals do not possess strong involutive bases. In 
particular, there is no point in invoking Proposition l5.7l for obtaining a strong basis. 
While we may surely obtain by elimination a subset H' such that Ic^ H' is a 
strong involutive basis of (le^ H), in general {H') C T. 

Example 9.10 Consider in the Weyl algebra W2 = li[x, y, dx, dy] the left ideal 
generated by the set = {l + x + y,dy — dx}- We take the semigroup order 
induced by the weight vector (—1,-1,1,1) and refined by a term order for which 
dy )^ dx ^ y >- X. Then the underlined terms are the leading ones. One easily 
checks that is a Grobner basis for this order. Furthermore, all variables are mul- 
tiplicative for each generator with respect to the Pommaret division and thus IF is 
a weak Pommaret basis, too. 

Obviously, the set F is neither a reduced Grobner basis nor a strong Pommaret 
basis, as 1 is a (multiplicative) divisor of dy. However, it is easy to see that the left 
ideal I — {T) does not possess a reduced Grobner basis or a strong Pommaret 
basis. Indeed, we have le^ I = Nq and thus such a basis had to consist of only a 
single generator; but I is not a principal ideal. <\ 

A special situation arises for the Janet division. Recall from Remark l276l that 
any finite set M C Ng is automatically involutively autoreduced with respect to 
the Janet division. Thus any weak Janet basis is a strong basis, if all generators 
have different leading exponents. If we follow the above outlined strategy of ap- 
plying Algorithm |3] to a homogenised basis and then dehomogenising the result, 
we cannot generally expect this condition to be satisfied. However, with a minor 
modification of the algorithm we can achieve this goal. 

Theorem 9.11 Let (T',*, -<) be a polynomial algebra of solvable type where -< 
is an arbitrary semigroup order Then every left ideal I (- V possesses a strong 
Janet basis for -<. 

Proof Assume that at some intermediate stage of Algorithm [3] the basis 7i con- 
tains two polynomials / and g such that le^^ (g) = Ic^^ (/) + lo, i. e. the leading 
exponents differ only in the first entry. If g — xqJ, we will find f = g after the 
dehomogenisation and no obstruction to a strong basis appears. Otherwise we note 
that, by definition of the lifted Janet division Jh, the homogenisation variable xq 
is non-multiplicative for /. Thus at some later stage the algorithm must consider 
the non-multiplicative product xqJ (if it was already treated, H would not be in- 
volutively head autoreduced). 

In the usual algorithm, we then determine the involutive normal form of the 
polynomial xof; the first step of this computation is to replace xof by xqJ — g. 
Alternatively, we may proceed instead as follows. The polynomial g is removed 
from the basis Ti, and replaced by x^f . Then we continue by analysing the invo- 
lutive normal form of g with respect to the new basis. Note that this modification 
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concerns only the situation that a multiplication by xq has been performed and 
that the basis Ti contains already an element with the same leading exponent as 
the obtained polynomial. 

If the final output Ti of the thus modified completion algorithm contains two 
polynomials / and g such that le^^ (g) and le_<^ (/) differ only in the first entry, 
then either g = Xq/ or / = x^g for some fc G N. Thus the dehomogenisation 
yields a basis 7i where all elements possess different leading exponents and TC 
is a strong Janet basis. Looking at the proof of Theorem 17.41 it is easy to see 
that this modification does not affect the correctness and the termination of the 
algorithm. As the Janet division is Noetherian, these considerations prove together 
with Proposition |9.4| the assertion. □ 

Note that our modification only achieves its goal, if we really restrict in Algo- 
rithm[3]to head reductions. Otherwise some other terms than the leading term in 
Xq / might be reducible but not the corresponding terms in /. Then we could still 
find after the dehomogenisation two generators with the same leading exponent. 

Example 9.12 Let us consider in the Weyl algebra W3 with the three variables x, 
y, z the left ideal generated by the set T — {dz — ydx, dy}. If we apply the usual 
involutive completion Algorithm[3](to the homogenisation JTC*) ), we obtain for the 
weight vector (—1, 0, 0, 1, 0, 0) refined by the degree reverse lexicographic order 
and the Janet division the following weak basis with nine generators: 

"Hi = { dj;, dy, dz, dxdz, dydz, ydx, yd^ + dz, yd^dz, yd^dz + 9f } . (28) 

As one easily sees from the last four generators, it is not a strong basis. 

Applying the modified algorithm for the Janet division yields the following 
basis with only seven generators: 

= { i9x + dydz, dy, dz, dxdz, dydz, ydx + dz, ydxdz + 9f } . (29) 

Obviously, we now have a strong basis, as all leading exponents are different. 

This example also demonstrates the profound effect of the homogenisation. 
A strong Janet or Pommaret basis of (JT) is simply given by = {dx, dy, dz} 
which is simultaneously a reduced Grobner basis. In (JTC*)) many reductions are 
not possible because the terms contain different powers of xq. However, this is a 
general problem of all approaches to Grobner bases for semigroup orders using 
homogenisation and not specific for the involutive approach. 

In this particular case, one could have applied the involutive completion algo- 
rithm directly to the original set T and it would have terminated with the minimal 
basis Ji, although we are using a order which is not a monoid order Unfortunately, 
it is not clear how to predict when infinite reduction chains appear in normal form 
computations with respect to such orders, so that one does not know in advance 
whether one may dispense with the homogenisation. < 
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10 Involutive Bases for Semigroup Orders II: Mora's Normal Form 

One computational disadvantage of the approach outHned in the previous section 
is that the basis Ti. in the homogenised algebra V is often much larger than the 
final basis Ti. in the original algebra V, as upon dehomogenisation generators may 
become identical. Furthermore, we have seen that it is difficult to prove the con- 
structivity of the lifted division Lh which limits the applicability of this technique. 
Finally, for most divisions we are not able to determine strong bases. 

An alternative approach for Grobner bases computations in the ordinary poly- 
nomial ring was proposed first by Greuel and Pfister ||29l and later independently 
by Grabe II27I I281: extensive textbook discussions are contained in fTS' Chapt. 4] 
and ll30l Sect. 1.6]. It allows us to dispense completely with computing in the ho- 
mogenised algebra V. Two ideas are the core of this approach: we modify the 
normal form algorithm using ideas developed by Mora for the computation of 
tangent cones and we work over a ring of fractions of V. We will now show that a 
generalisation to arbitrary polynomial algebras of solvable type and to involutive 
normal forms is possible and removes all the mentioned problems. 

The central problem in working with semigroup orders is that they are no 
longer well-orders and hence normal form computations in the classical form do 
not necessarily terminate. Mora |49 | introduced the notion of the ecart of a poly- 
nomial / as the difference between the lowest and the highest degree of a term in / 
and based a new normal form algorithm on it which always terminates. The main 
differences between it and the usual algorithm lie in the possibility to reduce also 
with respect to intermediate results (see Line 191 in Algorithm|5]below) and that it 
computes only a "weak" normal form (cf. Proposition [ToT|below). 

Mora's approach is valid only for tangent cone orders where the leading term is 
always of minimal degree. Greuel and Pfister [29 1 noticed that a slight modification 
of the definition of the ecart allows us to use it for arbitrary semigroup orders. So 
we set for any polynomial f €V \ {0} and any semigroup order -< 

ecart / = deg / - deg It^ / . (30) 

The extension of the Mora normal form to an involutive normal form faces 
one problem. As akeady mentioned, one allows here also reductions with respect 
to some intermediate results and thus one must decide on the assignment of mul- 
tiplicative variables to these. However, it immediately follows from the proof of 
the correctness of the Mora algorithm how this assignment must be done in or- 
der to obtain in the end an involutive standard representation with respect to the 
set Q (one should stress that this assignment is not performed according to some 
involutive division in the sense of Definition 12. 11 1. 

In Algorithm|5]below we use the following approach. To each member g of the 
set Q with respect to which we reduce we assign a set N[g\ of multiplicative in- 
dices. We write le^ g \n le^ h, if the multi index le^ h lies in the restricted cone of 
le^ g defined by N[g\. The set <S collects all generators g £ Q which have already 
been used for reductions and the set Af is the intersection of the corresponding 
sets of multiplicative indices. If a new polynomial h is added to Q, it is assigned 
as multiplicative indices the current value of M. 
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Algorithm 5 Involutive Mora normal form for a semigroup order ^ onV 
Require: polynomial / G P, finite set Q CV, involutive division L 
Ensure: involutive Mora normal form ft of / with respect to Q 

I: f- g^g 

2: for aWgeg do 

3: N[g]^ NL,ic^g{\e^g) 

4: end for 

5: AT ^ {!,..., n}; 5^0 

6: wiiile (/i / 0) A (3 g G (J : le^ 5 |jv le^ h) do 

7: choose g with ecart (; minimal among all g G g such that le^ g |jv le^ h 

8: if {g eg)A{g i S) tiien 

9: 5^5u{g}; ^ n iVfg] 

10: end if 

11: if ecart g > ecart h tlien 

12: g^^U{ft}; Aff/il^A/" 

13: end if 

14: /i le^ /i - le^ g; h~ i^^\t!^^g^, x'' * g 

15: end wliile 

16: return h 



Proposition 10.1 Algorithm\5\always terminates. Let {V — k[X], -<) be a poly- 
nomial algebra of solvable type (for an arbitrary semigroup order ~<) such that 
ii[X'] is a subring of V for any subset X' C X. Then the output h is a weak 
involutive normal form of the input f with respect to the set Q in the sense that 
there exists a polynomial u £ V with le^ u — such that the difference u-k f — h 
possesses an involutive standard representation 



and none of the leading exponents le^ g involutively divides le^ h. If ~< is a 
monoid order, then u = 1 and h is an involutive normal form in the usual sense. 

Proof As the proof is almost identical with the one for the non-involutive version 
of the Mora normal form given by Greuel and Pfister [29"301, we only sketch the 
required modifications; full details are given in [55J. For the termination proof no 
modifications are needed. For the existence of the involutive standard representa- 
tion one uses the same induction as in the non-involutive case and keeps track of 
the multiplicative variables. The key point is that if a reduction with respect to a 
polynomial .g G ^ \ 5 is performed, then this polynomial is multiplied only with 
terms which are multiplicative for a[\ g £ Q appearing in g. This ensures that in 
the end indeed each non-zero coefficient Pg is contained in ]k;[Xi g.^ (5)]. □ 

Remark 10.2 The assumption about V in Proposition 110. II is necessary, because 
the coefficients Pg in dsTl ) are the result of multiplications. While the above con- 
siderations ensure that each factor lies in ]k;[XL, 5,^(5)], it is unclear in a general 
polynomial algebra whether this remains true for their product. Simple examples 
for polynomial algebras of solvable type satisfying the made assumption are rings 




(31) 



see 
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of linear difference or differential operators. In the case of the Pommaret division, 
the assumption can be weaken a bit and every iterated polynomial algebra of solv- 
able type in the sense of Definition |4T| is permitted, too. <] 

We move now to a larger ring of fractions where all polynomials with leading 
exponent are units. In such a ring it really makes sense to call h a (weak) normal 
form of /, as we multiply / only by a unit. 

Proposition 10.3 Let (7^,-^, ~<) be a polynomial algebra of solvable type where ~< 
is a semigroup order Then the subset 

S^ = {f ^V\\e^f = Q} . (32) 

is multiplicatively closed and the left localisation = -k V is a well defined 
ring of left fractions. 

Proof Obviously, 1 G 5^ . If 1 + / and 1 + .g are two elements in S-^ , then the 
compatibility of the order ^ with the multiplication ★ ensures that their product 
is of the form (1 + /) ★ (1 + g) = 1 + /i with le^ h ^ Q. Hence the set S-^ is 
multiplicatively closed. 

As polynomial algebras of solvable type do not possess zero divisors, a suffi- 
cient condition for the existence of the ring of left fractions ★ P is that for all 
/ S and g the intersection {P*f)r\ {S^ -kg) is not empty (TS". Sect. 12. 1]. 
But this can be shown using minor modifications of our proof of Proposition l3.8l 

We first choose coefficients tq, sq E TZ such that in hi ~ r^g * / — sq/ * g 
the leading terms cancel, i. e. we have Ic^ hi -< Ic^ / + Ic^ g = Ic^ g. Then 
we compute with (the non-involutive form of) Algorithm |5] a weak normal form 
hi of hi with respect to the set = {f,g}- By Proposition II 0. II this yields a 
standard representation ui-khi — hi — (/>o * / + V-'o * .9 where le^ ui = 0. Assume 
that le^ -00 ^ 0. Then we arrive at the contradiction le^ {tpo -k g) h le^ g >- 
le^ hi = le^ (ui -khi). Thus le^ ^ 0. If hi = 0, then {ui ★ r^g — (j)o) -k f = 
{ui ★ So/ + V'o) * 5 and by the considerations above on the leading exponents 
ui k So/ + i>o eS^ so that indeed {V ^ f) Ci (S^ -kg) 

If hi ^ 0, we proceed as in the proof of Proposition [33] We introduce J^i = 
Tq U {hi\ and choose ri, si G such that in /i2 = rihi -k f — sif ★ hi the 
leading terms cancel. If we compute a weak Mora normal form /12 of /i2, then 
we obtain a standard representation U2 -kh2 — h2 = cj)i -k f + ^i -k g + pi -k hi 
where again le^ U2 = 0. The properties of a standard representation imply now 
that le^ -01 + le^ g d: le^ h2 and le^ pi + le^ hi ^ le^ /i2- Together with the 
inequalities le^ /12 -< le^ / + le^ hi = le^ hi -< le^ g this entails that both 
le^ -01 ^ and le^ pi -< 0. Thus for /12 = we have found (f> E V and -0 € 
such that (j>-kf = tpkg.lfh2 7^ 0, similar inequahties in the subsequent iterations 
ensure that we always have S 5^ . □ 

As any localisation of a Noetherian ring is again Noetherian, is Noetherian, 
if V is so. One sees immediately that the units in are all those fractions where 
not only the denominator but also the numerator is contained in . Given an 
ideal 2 C "P^, we may always assume without loss of generality that its generated 
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by a set c P of polynomials, as multiplication of a generator by a unit does 
not change the span. Hence in all computations we will exclusively work with 
polynomials and not with fractions. 

As all elements of are units in V^, we may extend the notions of leading 
term, monomial or exponent: if / G V^, then we can choose a unit u e 5^ with 
Ic^ M = 1 such that ut^t / e P is a polynomial; now we define le^ / — le^ (u ★ /) 
etc. One easily verifies that this definition is independent of the choice of u. 

Following Greuel and Pfister |30|, one can now construct a complete theory 
of involutive bases over . Definition 19.21 of Grobner and involutive bases can 
be extended without changes from the ring P to V^. Theorem 14.71 on the exis- 
tence of Grobner bases generalises to P^, as its proof is only based on the leading 
exponents and a simple normal form argument remaining valid due to our consid- 
erations above. 

Note that even if the set Q is involutively head autoreduced, we cannot con- 
clude in analogy to Proposition 15.131 that the involutive Mora normal form is 
unique, as we only consider the leading term in Algorithm |5] and hence the lower 
terms in h may still be involutively divisible by the leading term of some generator 
g € G- However, Theorem l5.4l remains valid. 

Theorem 10.4 Let {P — k[X],*, -<) be a polynomial algebra of solvable type 
(for an arbitrary semigroup order such that k[X'] is a subring of P for any 
subset X' C X. Furthermore, let L be a constructive Noetherian division. For a 
finite set T (ZP of polynomials let X = {J-) be the left ideal generated by it in the 
localisation P^. If we apply Algorithm\3\with the involutive Mora normal form 
instead of the usual one to the set T, then it terminates with an involutive basis of 
the ideal X. 

Proof The termination of Algorithm |3] under the made assumptions was shown 
in Proposition 17.21 and Theorem 17.41 One easily verifies that their proofs are not 
affected by the substitution of the normal form algorithm, as they rely mainly 
on Theorem 15.41 and on the fact that the leading term of the normal form is not 
involutively divisible by the leading term of any generator Both properties remain 
valid for the Mora normal form. □ 

Remark 10.5 Note that Theorem ll0.4| guarantees the existence of strong involutive 
bases. Due to the extension to P^, Example 19.101 is no longer a valid counter 
example. As the first generator in T is now a unit, we find that (JT) = P^ and {1} 
is a trivial strong Pommaret basis. <\ 

Example 10.6 We continue Example l9.12l Following the approach given by The- 
orem ll0.4l we immediately compute as Janet basis of {T) (over P^) the minimal 
basis 7^3 ~ {dx, dy, dz}. Obviously, it is considerably smaller than the bases ob- 
tained with Lazard's approach (over P). This effect becomes even more profound, 
if we look at the sizes of the bases in the homogenised Weyl algebra; both Hi and 
7^2 consist of 21 generators. < 
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11 Involutive Bases over Rings 

Finally, we consider the general case that V — Ti-lxi, . . . polynomial 
algebra of solvable type over a (left) Noetherian ring TZ. In the commutative case, 
Grobner bases for such algebras have been studied in [25','59l (see fT, Chapt. 4] for 
a more extensive textbook discussion); for PBW extensions (recall Example l3.5l ) a 
theory of Grobner bases was recently developed in 1261 . We will follow the basic 
ideas developed in these references and assume in the sequel that linear equations 
are solvable in the coefficient ring TZ, meaning that the following two operations 
can be effectively performed: 

(i) given elements s, ri, . . . , rj; G TZ, we can decide whether s G {ri, ■ ■ ■ ,rk)'R. 
(the left ideal in TZ generated by ri , . . . , r^,); 

(ii) given elements ri , . . . , ^ TZ, we can construct a finite basis of the module 

Syz(ri, . . . , Tfe) of left syzygies siri H h 5^.?'^ = 0. 

The first operation is obviously necessary for the algorithmic reduction of poly- 
nomials with respect to a set c V. The necessity of the second operation will 
become evident later Compared with the commutative case, reduction is a more 
complicated process, in particular due to the possibility that in the commutation 
relations ^ for the multiplication in V the maps may be different from the 
identity on TZ and the coefficients r^^ unequal one. 

Let Q C T-" be a finite set. We introduce for any polynomial f ^ P the sets 
Gf ^ {geg\ le^ 5 I le^ /} and 

Gf = {xf^^glgeGf A^i^le^f~\e^gA le^ (x'' * g) = Ic^ /} (33) 

Note that the last condition in the definition of Q f is redundant only, if the coef- 
ficient ring TZ is an integral domain. Otherwise it may happen that \Qf\ < \Gf\, 
namely if p^(r)r^^ = where Im^ g = rx'^ . The polynomial / is head reducible 
with respect to Q, if Ic^ g G (Ic^ (note that we use Qf here so that the re- 
duction comes only from the leading terms and is not due to some zero divisors as 
leading coefficients). Involutive head reducibility is defined analogously via sets 
Q fx and Q fx where only involutive divisors with respect to the division L on Ng 
are taken into account, i. e. we set 

Gf,L = {5 e g I le^ / G Clm, e(le^ g)} . (34) 

Thus the set Q is involutively head autoreduced, if Ic^ g ^ (Ic^ {Gg,L \ {g}))TZ 
for all polynomials g ^ G- This is now a much weaker notion than before; in 
particular, Lemma l5.12l is no longer valid. 

Definition 11.1 Let X C T' be a left ideal in the polynomial algebra (T',*, -<) 
of solvable type over a ring TZ in which linear equations can be solved. A finite 
set G T' is a Grobner basis of X, if for every polynomial / G I the condition 
Ic^ / G (Ic^ G})ti is satisfied. The set G is a weak involutive basis/or the involu- 
tive division L, if for every polynomial / G I the condition Ic^ / G (Ic^ Gf,L)Ti 
is satisfied. A weak involutive basis is a strong involutive basis, if every set Gf,L 
contains precisely one element. 
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It is easy to see that the characterisation of (weak) involutive bases via the exis- 
tence of involutive standard representations (Theorem 15 .4b remains valid. Indeed, 
only the first part of the proof requires a minor change: the polynomial /i is now 
of the form /i = / — J^hei-if l ^'^^ where the coefficients E TZ sue chosen 
such that le^ /i -< le^ /. 

Clearly, a necessary condition for the existence of Grobner and thus of (weak) 
involutive bases for arbitrary left ideals T C V is that the algebra is a (left) 
Noetherian ring. As we have seen in Section ID this assumption becomes non- 
trivial, if the coefficient ring TZ is not a field. In this section, we will assume 
throughout that T' is a polynomial algebra of solvable type over a left Noetherian 
ring TZ with centred commutation relations (cf . Definition l4.3l l so that Theorem l4.4l 
asserts that V is left Noetherian, tooH A very useful side effect of this assumption 
is that the scalars appearing in the commutation relations ^ are units and thus not 
zero divisors which is important for some arguments. 

Example 11.2 As in the previous two sections, we cannot generally expect strong 
involutive bases to exist. As a simple concrete example, also demonstrating the 
need of the second assumption on TZ, we consider in k[a;, y\[z\ (with the ordinary 
multiplication) the ideal I generated by the setjr = {x^z — l,y^z+l}. Obviously, 
both generators have the same leading exponent [1]; nevertheless none is reducible 
by the other one due to the relative primeness of the coefficients. Furthermore, the 
syzygy S — x^e.i — y^ei G W^x. y]"^ connecting the leading coefficients leads to 
the polynomial + € I. It is easy to see that a Grobner and weak Janet basis 
of I is obtained by adding it to T. A strong Janet basis does not exist, as none of 
these generators may be removed from the basis. <\ 

This example shows that simply applying our completion Algorithm [3] will 
generally not suffice. Obviously, with respect to the Janet division z is multiplica- 
tive for both elements of T so that no non-multiplicative variables exist and thus 
it is not possible to generate the missing generator by multiplication with a non- 
multiplicative variable. We must substitute in Algorithm |3] the involutive head au- 
toreduction by a more comprehensive operation^ 

Definition 11.3 Let J- C T' be a finite set and L an involutive division. We con- 
sider for each f E T the syzygies ^ f^fj ^ 1*^^ / — connecting the leading 

* The case of an iterated polynomial algebra of solvable type (cf. Definition 14. It will be 
considered in Part II, after we have developed a syzygy theory for involutive bases. 

' In the classical case of commutative variables over a coefficient field, it is not difficult to 
show that for any finite set the syzygy module Syz(lm^ T) of the leading monomials can 
be spanned by binomial generators corresponding to the S-polynomials in the Buchberger 
algorithm. In Part II we will show that in any such syzygy at least one component contains 
a non-multiplicative variable, so that implicitly the involutive completion algorithm also 
runs over a generating set of this syzygy module. When we move on to coefficient rings, 
it is well-known that additional, more complicated syzygies coming from the coefficients 
must be considered. For these we can no longer assume that one component contains a non- 
multiplicative variable. Hence partially we must follow the same approach as in the gen- 
eralisation of the Buchberger algorithm and this leads to the notion of 7?,-saturation where 
some syzygies not reachable via non-multiplicative variables are explicitly considered. 
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coefficients of the elements of the set Tfx- The set J- is involutively 7?.- saturated 
for the division L, if for any such syzygy S the polynomial ^ f£fj z, ^ff possesses 
an involutive standard representation with respect to J-. 

For checking involutive 7?,-saturation, it obviously suffices to consider a finite 
basis of each of the finitely many syzygy modules Syz(lc^ -^/.l) so that such 
a check can easily be performed effectively. An element f G !F is involutively 
head reducible by the other elements of JF, if and only if Syz(lc^ ^f,L) contains 
a syzygy with Sf — 1. For this reason it is easy to combine an involutive TZ- 
saturation with an involutive head autoreduction leading to Algorithm|6] 

The for loop in Lines /5-13/ takes care of the involutive head autoreduction 
(the call HeadReducei.^(/, 7i) involutively head reduces / with respect to the 
set H \ {/} but with multiplicative variables determined with respect to the full 
set Ti — cf. Remark |S!9] l. The for loop in Lines 111-221 checks the involutive TZ- 
saturation. Each iteration of the outer while loop analyses from the remaining 
polynomials (collected in S) those with the highest leading exponent. The set S is 
reset to the full basis, whenever a new element has been put into TL; this ensures 
that all new reduction possibilities are taken into account. In Line /15/ it does not 
matter which element / G 5^ we choose, as the set Ti'j ^ depends only on Ic^ / 
and all elements of Si, possess by construction the same leading exponent v. 

Proposition 11.4 Under the made assumptions about the polynomial algebra V, 
Algorithm^terminates for any finite input set J- (Z V with an involutively TZ- 
saturated and head autoreduced set TL such that {H) = {J-). 

Proof The correctness of the algorithm is trivial. The termination follows from the 
fact that both TZ and INq are Noetherian. Whenever we add a new polynomial h 
to the set H', we have either that le^ h ^ (le^ or Ic^ h ^ (Ic^ H'l^ j^)ti. 

As neither in Ng nor in TZ infinite ascending chains of ideals are possible, the 
algorithm must terminate after a finite number of steps. □ 

An obvious idea is now to substitute in the completion Algorithm[3]the involu- 
tive head autoreduction by an involutive 7?.- saturation. Recall that Proposition l7.2l 
(and Corollarv l7.3b was the crucial step for proving the correctness of Algorithmic] 
Our next goal is thus to show that under the made assumptions for involutively TZ- 
saturated sets local involution implies weak involution. 

Proposition 11.5 Under the made assumptions about the polynomial algebra V, 
a finite, involutively TZ-saturated set J- <Z V is weakly involutive, if and only if it 
is locally involutive. 

Proof We first note that Proposition l7.2l remains true under the made assumptions. 
Its proof only requires a few trivial modifications, as all appearing coefficients 
(for example, when we rewrite x^^ — > x^^~^^ * Xj) are units in the case of centred 
commutation relations and thus we may proceed as for a field. Hence if T is locally 
involutive, then T = {T) = (JT) implying that any polynomial g E 2 may be 
written in the form g — J^fej^Pf * / ^i'^^ Pf G 7?.[Xl,jc-.^(/)]. Furthermore, 
it follows from this proof that for centred commutation relations we may assume 
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Algorithm 6 Involutive 7?.-saturation (and head autoreduction) 

Require: finite set T <ZV, involutive division L on Nq 

Ensure: Involutively 7^-saturated and head autoreduced set H. with (7i) = {T) 

\: T\ S ^ T 

2: while S / do 

3: 1/ «— max^ le^ S\ S,j ^ {f £ H \ le^ f ~ ^} 

4: 5 ^ 5 \ S,/, H' 

5: for all / G 5^, do 

6: ft ^ HeadReducei,^^ (/, 7i) 

7: if ///i then 

8: 5. ^5,A{/}; n'^n'Mf} 

9: if /i 7^ then 

10: H' ^H'U {h} 

11: end if 

12: end if 

13: end for 

14: if 5^/0 then 

15: choose f £ and determine the set 

16: compute basis S of Syz(lc^ 

17: for all S = E/gw'. ^f^f G ^ do 

18: /l ^ NormalFormL,-<(E/eH' Sff,H') 

19: if /i 7^ then 

20: H' ^H'U {h} 

21: end if 

22: end for 

23: end if 

24: if-WVWthen 

25: H^H'; S^H 

26: end if 

27: end while 

28: return H 



that the polynomials P/ satisfy le^ (P/ * f) = Ic^ P/ + le^ /. We are done, if we 
can show that they can be chosen such that additionally Ic^ (Py * f ) Ic^ g, i. e. 
such that we obtain an involutive standard representation of g. 

If the representation coming out of the proof of Proposition l7.2l alreadv satisfies 
this condition on the leading exponents, nothing has to be done. Otherwise we set 

i/ = max^ {le^ (P/ /) | / e T} and ^ {f e T \ le^ (P/ ★ /) = ly}. As by 
construction i/ € ClfeJ^ ^i.io^ /), the properties of an involutive division 

imply that we can write T„ ^ {fi, . . . , fk} with le^ /i | le^ /2 | • • ■ | le^ fk and 
hence JT^ C J-f^.L- Since we have assumed that le^ (P/ -*:/) — le^ P/ + le^ /, 
we even find C J-j^.l- 

By construction, the equaUty '^ffzjr Ic^ (P/ * /) = holds. If we now set 
Im^ / = rfX^f and Im^ P/ — SfX^f , then we obtain under the made assump- 
tions: Ic^ (P/ ★ /) = SfPf,j{rf)rf,j^j = [sfPnj(rf)r^ji,j\rf and hence the 
above equality corresponds to a syzygy of the set Ic^ ^fk,L- As the set T is invo- 
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lutively 7?.-saturated, there exists an involutive standard representation 

fe 

withQ/ e k[XL,^,^(/)] and le^ (Q/*/) = le^ Q/ + le^ / -< Vf^. 

Introducing now the polynomials Q'^ = Qf — [s fP^f{rf)rf^^,j^.'\x'^^k~^f for 
f E and = Q/ otherwise, we get the syzygy J^fej^ Q'f * / = 0. If we set 
P'f ^ Pf - cj^x^-^ff^ *Q'j: with Cf = Pv-u;^ {sfPt^firf)r,,f^f)p^firf)r^^^^, 
then, by construction, g — J^feJ^P'f * f another involutive representation of 
the polynomial g with i^' = max^ {le^ (Pj- * f ) \ f E J-} ^ ly. 

Repeating this procedure for a finite number of times obviously yields an invo- 
lutive standard representation of the polynomial g. As g was an arbitrary element 
of the ideal T = {T), this implies that T is indeed weakly involutive. □ 

Theorem 11.6 Let V be a polynomial algebra of solvable type satisfying the made 
assumptions. If the subalgorithm InvHeadAutoReduce/,^^ is substituted inAlgo- 
rithm\3\by Algorithm^ then the completion will terminate with a weak involutive 
basis of I = {J-) for any finite input set !F G V such that the monoid ideal le^ I 
possesses a weak involutive basis. 

Proof The correctness of the modified algorithm follows immediately from Propo- 
sition [TT3] For the termination we may use the same argument as in the proof of 
Theorem l7.4l as it depends only on the leading exponents. □ 

12 Conclusions 

We studied involutive bases for a rather general class of non-commutative polyno- 
mial algebras. Our approach was very closely modelled on that of Kandry-Rody 
and Weispfenning |41| and subsequently Kredel [42|. However, we believe that 
the third condition in Definition [TT] (the compatibility between term order -< and 
non-commutative product *) is more natural than the stricter axioms in [411. In 
particular, we could not see where Kandry-Rody and Weispfenning needed their 
stricter conditions, as all their main results hold in our more general situation, as 
later shown by Kredel. 

Comparing with ||2l [T3ll42l|45l . one must say that the there used approach is 
more constructive than ours. More precisely, all these authors specify the non- 
commutative product via commutation relations and thus have automatically a 
concrete algorithm for evaluating any product. As we have seen in the proof of 
Proposition l3.3l the same data suffices to fix our axiomatically described product, 
but it does not provide us with an algorithm. However, we showed that we can 
always map to their approach via a basis transformation. 

We showed that the polynomial algebras of solvable type form a natural frame- 
work for involutive bases. This is not surprising, if one takes into account that the 
main part of the involutive theory happens in the monoid Ng and the decisive third 
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condition in Definition 13. II of a polynomial algebra of solvable type ensures that 
its product * does not interfere with the leading exponents. 

We extended the theory of involutive bases to semigroup orders and to poly- 
nomials over coefficient rings. It turned out that the novel concept of a weak in- 
volutive basis is crucial for such generalisations, as in both cases strong bases 
rarely exist. These weak bases are still Grobner bases and involutive standard rep- 
resentations still exist (though they are no longer unique). It seems that in such 
computations the Janet division has a distinguished position, as by Theorem 19. Ill 
strong Janet bases always exist. If one is only interested in using Algorithm|3]as an 
alternative to Buchberger's algorithm, weak bases are sufficient. However, most of 
the more advanced applications of involutive bases studied in Part II will require 
strong involutive bases. 

Concerning involutive bases over rings, we will study in Part II the special 
case that the coefficient ring is again a polynomial algebra of solvable type. Using 
the syzygy theory that will be developed there, we will be able to obtain stronger 
results and a "purely involutive" completion algorithm. The current approach con- 
tains hidden in the concept of 7?,-saturation parts of the Buchberger algorithm for 
the construction of Grobner bases over rings. 

Definition |2T|represents the currently mainly used definition of an involutive 
division. While it appears quite natural, one problem is that in some sense too many 
involutive divisions exist, in particular rather weird ones with unpleasant proper- 
ties. This effect has lead us to the introduction of such technical concepts like 
continuity and constructivity. One could imagine that there should exist a stricter 
definition of involutive divisions that automatically ensures that Algorithmic ter- 
minates without having to resort to these technicalities. 

Most of these weird divisions are globally defined and multiplicative indices 
are assigned only to finitely many multi indices. Such divisions are obviously of 
no interest, as more or less no monoid ideal possesses an involutive basis for them. 
One way to eliminate these divisions would be to require that for every q G No 
the monoid ideal (No)>g = {ly € Nq | q < has an involutive basis. All the 
involutive divisions used in practice satisfy this condition, but it is still a long way 
from this simple condition to the termination of Algorithm|2] 

We did not discuss the efficiency of the here presented algorithms. Much of 
the literature on involutive bases is concerned with their use as an alternative ap- 
proach to the construction of Grobner bases. In particular, recent experiments by 
Gerdt et al. ||23 | comparing a specialised C/C++ program for the construction of 
Janet bases with the Grobner bases package of SINGULAR |[3T| indicate that the 
involutive approach is highly competitive. This is quite remarkable, if one takes 
into account that SINGULAR is based on the results of many years of intensive 
research on Grobner bases by many groups, whereas involutive bases are still very 
young and only a few researchers have actively worked on them. The results in 
Part II will offer some heuristic explanations for this observation. 

Finally, we mention that most of the algorithms discussed in this article have 
been implemented (for general polynomial algebras of solvable type) by M. Haus- 
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dorf 13311341 in the computer algebra system MuPAD^ The implementation does 
not use the simple completion Algorithm[3]but a more optimised version yielding 
minimal bases developed by Gerdt and Blinkov [21] , It also includes the modified 
algorithm for determining strong Janet bases in local rings. 

A Term Orders 

We use in this article non-standard definitions of some basic term orders. More 
precisely, we revert the order of the variables: our definitions become the standard 
ones, if one transforms (xi, . . . , Xn) {xn, ■ ■ ■ , xi). The reason for this reversal 
is that this way the definitions fit better to the conventions in the theory of invo- 
lutive systems of differential equations. Furthermore, they appear more natural in 
some applications like the determination of the depth in Part II. 

A term order -< is for us a total order on the set T of all terms x^ satisfying 
the following two conditions: (i) 1 < t for all terms t £ T and (ii) s ^ t implies 
r ■ s ^ r -t for all terms r,s,t £ T. If a term order fulfils in addition the condition 
that s ^ t whenever deg s < deg t, it is called degree compatible. As T and Ng 
are isomorphic as monoids, we may also speak of term orders on Nq . In fact, most 
term orders are defined via multi indices. 

A more appropriate name for term orders might be monoid orders, as the two 
conditions above say nothing but that these orders respect the monoid structure 
of T. A more general class of (total) orders are semigroup orders where we skip 
the first condition, i. e. we only take the semigroup structure of T into account. It 
is a well-known property of such orders that they are no longer well-orders. This 
implies in particular the existence of infinite descending sequences so that normal 
form algorithms do not necessarily terminate. 

The lexicographic order is defined by x'^ ^lex x'^, if the last non-vanishing 
entry of /i — is negative. Thus X2X3 -(lex 3:1X3. With respect to the reverse 
lexicographic order x'^ ^leviex x'^, if the first non-vanishing entry of /i — 1/ is 
positive. Now we have xix^ -<reviex 2:2X3. However, -(reviex is only a semigroup 
order, as it violates the first condition: xi ^leviex 1- Degree compatible versions of 
these orders exist, too. x^ -<degiex x'^, if < or if = and x^ -<iex x". 
Similarly, x^ -<degi-eviex x"^, if < or if = andx^ -<reviex x". Obviously 
^degreviex IS a term Order It possesses the following useful characterisation which 
is easy to prove. 

Lemma A.l Let -< be a degree compatible term order such that the condition 
It^ / G (xi, . . . , Xfc) is equivalent to f £ (xi, . . . , Xk) for every homogeneous 
polynomial / £ "P. Then -< is the degree reverse lexicographic order ^deereviex- 

We say that a term order respects classes, if for multi indices /i, v of the same 
length els fi < chv implies x^ -< x'^. It is now easy to see that by Lemma [ATl on 
terms of the same degree any class respecting term order on T coincides with the 
degree reverse lexicographic order If we consider free polynomial modules, class 
respecting orders have the same relation to TOP lifts |[T| of ^degreviex- 
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